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Abstract 


An  instrumentation  package  was  developed  for  measuring  the  electrical 
conductivity  of  flowing  MHD  plasmas  with  four-pin  probes.  This  system  was 
then  used  to  investigate  the  sensitivity  of  the  indicated  conductivity  to 
the  flow  velocity  and  the  pin  orientation  in  a subsonic,  atmospheric- 
pressure  argon  flow  at  2000  K seeded  with  0.1%  potassium. 

Experimental  investigation  of  the  current-voltage  characteristics  of 
such  probes  demonstrated  the  linearity  of  the  charactistics  for  applied 
voltages  of  about  1 volt  across  the  outer  pins,  and  pointed  out  the  need 
for  a.c.  excitation  of  the  probes.  When  excited  at  audio  frequencies 
(typically  20  Khz)  well  above  the  highest  frequencies  of  the  conductivity 
fluctuations  in  the  flow,  the  final  instrumentation  system  was  found  to  be 
well  suited  for  time-resolved  relative  conductivity  fluctuation  measurements 
and  for  continuous  display  of  the  time-averaged  indicated  conductivity, 

'^indic. 

The  absolute  accuracy  of  conductivity  measurements  using  four-pin  probes 

in  practical  MHD  flows  appears  to  be  limited  by  relatively  large  "flow 

effects."  For  example,  a decrease  in  approximately  50%  was 

observed  as  the  flow  velocity  increased  from  10  m/sec  to  100  m/sec.  Also, 

o.  decreased  as  the  pins  were  rotated  from  the  "end-on"  to  a sidewise 

indic 

orientation. 
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Chapter  1 


‘ INTRODUCTION 

1 . 1 Motivation  for  the  Research 

The  work  described  in  this  report  stems  from  a broader  interest  in 
studying  how  the  performance  of  MHD  generators  is  affected  by  spatial  and 
temporal  nonuniformities  which  exist,  for  various  reasons  and  to  varying 
degrees,  in  all  real  MHD  plasmas.  Some  of  the  mechanisms  which  may  pro- 
duce nonuniformities  in  the  plasma  properties  (e.g.,  electrical  conducti- 
vity, temperature,  electron  density,  velocity)  are  plasma  instabilities, 
poor  mixing,  imperfect  combustion,  and  three-dimensional  fluid  mechanical 
flow  effects  such  as  turbulence  and  boundary  layers.  At  the  present  time, 
an  Important  goal  for  researchers  in  the  field  of  MHD  power  generation  is 
to  develop  accurate  mathematical  models  for  predicting  the  performance  of 
proposed  large-sized  generators.  As  such  models  grow  more  refined,  it 
will  be  necessary  that  they  take  into  account  the  presence  of  nonuniformi- 
ties in  some  way,  via  the  parameters  which  characterize  these  nonuniformi- 
ties (such  as,  for  example,  statistical  measures  including  mean  and  rms 
values  of  fluctuating  conductivity,  probability  density  functions,  charac- 
teristic frequencies  or  frequency  spectra  of  fluctuations  in  plasma  prop- 
erties, and  many  others).  Such  information  will  be  required  as  input 
Information  to  the  theoretical  models.  For  this  b;^aic  reason,  MHD  research- 

\ 

ers  will  benefit  by  being  able  to  measure  experimentally  the  plasma  s non- 
uniform  and  fluctuating  properties,  with  three-dimensional  spatial  resolution 
and  adequate  time  resolution,  in  laboratory-scale  experiments.  In  summary, 
there  is  a need  for  spatially  and  temporally  resolved  measurements  of  MHD 
plasma  properties  to  help  formulate  and  provide  quantitative  input  to  more 
accurate  performance-predicting  models. 

Of  particular  importance  in  this  regard  is  the  development  of  a diag- 
nostic technique  for  measuring  the  electrical  conductivity  (or,  somewhat 
equivalently  under  most  conditions,  the  electron  density).  The  desirability  i 

of  a method  having  three-dimensional  spatial  resolution  stems  from  the  i 

essentially  three-dimensional  nature  of  many  plasma  nonuniformities,  from  | 

the  near-total  lack  of  any  actual  published  three-dimensional  fluctuations  j 

i 

f 

! 
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data  (e.g.,  measured  values  of  conductivity  fluctuations  at,  or  in  a small 
region  surrounding,  a point  in  the  flow  field)  up  to  the  present  time,  as 
well  as  from  the  need  to  assess  the  validity  of  more  traditional  optical 
diagnostics  which  give  only  two-dimensional  resolution  and  fail  to  resolve 
nonuniformities  along  the  optical  path.  In  addition  to  spatial  resolution, 
time-resolved  measurement  techniques  are  valuable  because  they  permit  fre- 
quency spectral  analyses  of  fluctuations  to  be  made  (which  often  illuminate 
the  nature  and  origin  of  the  nonuniformities) , and  offer  the  possibility 
(when  coupled  with  suitable  data-processing  equipment)  of  generating  any 
sort  of  quantitative  statistical  information  which  might  be  required  by 
theoretical  models,  such  as  the  theory  of  Reference  [1]. 

1 . 2 Review  of  Possible  Diagnostic  Techniques 

Measuring  techniques  for  electrical  conductivity  or  electron  density 
which  fulfill  the  requirements  for  both  three-dimensional  spatial  resolu- 
tion and  temporal  resolution  appear,  at  the  present  time,  to  be  restricted 
to  probe  methods,  viz,  the  four-electrode  (conductive  type)  conductivity 
probe,  the  inductive  RF  probes  (several  different  designs  exist),  and 
Langmuir  probes.  These  probes  are  described  in  Reference  [2]  and  will  be 
discussed  briefly  here.  First,  however,  we  mention  that,  while  there  are 
standard  spectroscopic  techniques  (such  as  the  absolute  continuum  intensity 
method  and  the  absolute  line  intensity  method  for  N^)  which  do  provide 
suitable  time  response,  their  spatial  resolution  is  inherently  limited  to 
two  dimensions.  Also,  there  are  laser  scattering  methods  (in  particular, 
Thomson  scattering  from  free  electrons  for  determining  N^)  which  are 
capable  of  3-D  spatial  resolution;  such  techniques  are  only  poorly  devel- 
oped at  the  present  time,  however,  and  for  MHD  plasmas  the  time  resolution 
is  likely  to  be  inadequate,  due  to  the  very  low  scattering  cross  sections 
and  consequent  long  times  required  for  photon  counting. 

The  four-electrode  probe  (or  four-pin  probe)  is  the  subject  of  this 
report  and  will  be  thoroughly  described  later;  it  is  perhaps  the  most 
"direct"  technique,  in  that  it  is  based  on  actually  passing  current  through 
the  plasma,  obtaining  measurements  of  both  the  current  density  J and 
electric  field  E,  and  deducing  the  conductivity  o directly  from  Ohm's 
Law,  J = aE. 
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The  several  types  of  inductive  (RF)  probes  which  have  been  used  to 
measure  o [2,3,4]  rely  on  the  fact  that  the  impedance  (or  Q)  of  a 
small  RF  coil  immersed  in  a plasma,  or  the  trans- impedance  between  two  or 
more  such  colls,  changes  when  the  conductivity  of  the  surrounding  medium 
changes.  The  colls  are  usually  enclosed  in  an  insulating  ceramic  tube 
and  swept  through  a plasma  stream  quickly  enough  to  avoid  melting.  The 
effects  of  nearby  metallic  conductors  tend  to  "swamp  out"  the  effects  of 
the  plasma,  whose  conductivity  is  orders  of  magnitude  lower;  in  addition, 
to  interpret  the  data  obtained  with  these  Inductive  probes  requires  assump- 
tions about  the  fluid  mechanical  perturbations  caused  by  the  presence  of 
the  probe.  These  facts  tend  to  limit  the  accuracy  and  the  spatial  resolu- 
tion of  inductive  probes,  although  it  appears  difficult  to  characterize 
this  limitation  quantitatively. 

The  electron  concentration,  N , can  be  determined  under  MHD  condi- 

e 

tions  by  suitably  designed  and  calibrated  Langmuir  probes.  Since  a = 

2 — — 
e the  collision  frequency  can  be  calculated  fairly 

accurately  if  the  temperature  is  known  to  only  modest  accuracy,  the  Lang- 
muir probe  measurement  together  with  a rough  temperature  measurement 

will  yield  the  electrical  conductivity.  The  theoretical  and  experimental 
work  of  Clements  [5]  has  demonstrated  the  feasibility  of  measurements 

by  Langmuir  probes  operated  in  the  highly  negatively-biased  ion-collecting 
mode.  However,  the  theory  which  relates  the  calculated  to  the  mea- 

sured probe  bias  voltage  and  collected  ion  current  is  only  approximate  at 

the  present  time.  Thus,  in  order  to  obtain  accurate  N measurements,  one 

e 

will  be  forced  to  calibrate  the  Langmuir  probe's  response  by  immersing  the 
probe  in  plasmas  in  which  is  accurately  known  by  some  otner  means. 

Aside  from  this  difficulty  vis-a-vis  calibration,  the  Langmuir  probe  ap- 
pears to  offer  the  possibility  of  excellent  spatial  and  temporal  resolution. 
There  are  practical  problems  with  such  probes,  however,  among  which  are  the 
cooling  of  the  probe  and  the  deposition  of  seed  materials  or  other  sub- 
stances in  the  flow  onto  the  cold  probe. 
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1.3  The  Four-Electrode  Conductivity  Probe 

This  report  describes  experimental  work  directed  towards  assessing 
the  suitability  of  the  four-electrode  conductivity  probe  (or  "four-pin 
probe")  for  three-dimensional  spatially-resolved  and  time-resolved  elec- 
trical conductivity  measurements  in  MHD  plasma  flows.  In  this  section  we 
present  a brief  description  of  the  construction  and  operation  of  these 
probes . 

The  four-pin  probe  consists  of  four  exposed  metal  electrodes  in  the 
form  of  short  tips  of  bare  wire  emerging  from  two  parallel  double-bore 
ceramic  tubes,  as  shown  in  Fig.  1.1.  The  probe  tips  are  immersed  in  a 
conducting  medium  whose  conductivity,  a,  is  to  be  determined.  An  exter- 
nal power  supply  (such  as  a battery  or  a more  elaborate  oscillator)  causes 
a current  i to  flow  between  the  two  outer  electrodes,  and  the  magnitude 
of  this  current  is  measured,  for  example,  by  measuring  the  voltage 
appearing  across  a resistor  ^gjiunt  series  with  the  outer  electrodes. 
Simultaneously,  one  measures  the  potential  difference  ‘^^inner  ^P~ 

pears  between  the  two  inner  pins;  this  is  done  using  a very  high  impedance 
voltmeter  so  that  no  current  flows  into  either  of  these  electrodes.  If 
conditions  are  such  that  the  probe  characteristics  are  linear,  i.e.,  that 

the  current  i and  voltage  difference  6V.  are  directly  proportional 

° inner 

to  each  other,  we  may  write 


inner 


= (proportionality!  /resistivity  of|  ^ 

' constant  / \ the  medium  / 


inner 


k i (J) 

p ' 0 ' 


(1.1) 


The  theoretical  justification  for  this  relationship  is  simply  that  in  an 
ohmic  medium,  such  as  we  assume  the  plasma  to  be  in  the  absence  of  a mag- 
netic field,  the  current  density  (and  hence  total  current,  1)  is  propor- 
tional to  the  potential  difference  between  two  measuring  stations  (viz, 
the  two  inner  pins) . The  reason  for  using  the  inner  pins  (which  draw  no 
current)  rather  than  the  outer  pins  (which  do  draw  current)  for  this  p'iien- 
tial  difference  measurement  is  that  plasma  sheath  voltage  drops  occur  near 


rfkii*-.  


the  current-carrying  electrodes.  The  magnitudes  of  these  drops  are  of 
the  same  order  or  greater  than  the  ohmic  voltage  drop  through  the  plasma, 
and  furthermore  are  unequal  at  the  cathode  and  anode.  Thus  the  measured 
voltage  difference  between  the  outer  pins  is  not  proportional  to  the  true 
ohmic  potential  drop  across  the  plasma.  However,  at  the  floating  inner 
pins,  the  plasma  sheath  drops  will  be  smaller  and  presumably  equal  (or 
nearly  so)  at  each  pin,  so  that,  when  the  potential  difference  between 
the  inner  pins  is  measured,  the  sheath  drops  subtract  out,  leaving  only 
the  ohmic  contribution. 

The  schematic  portrayal,  in  Fig.  1.2,  of  the  potential  distribution 
along  a line  joining  the  four  pins  may  help  to  clarify  the  preceding  ex- 
planation. In  this  figure  AV^  p represents  the  potential  drop  immedi- 
ately outside  the  metal  electrode  surface  due  to  the  fact  that  there  is 
some  work  function  potential  which  electrons  must  overcome  in  order  to 
leave  the  metal.  Also,  ^^sheath  represents  the  plasma  sheath  voltage 
drop,  and,  as  mentioned  previously,  there  is  an  ohmic  potential  drop  of 

k pi,  or  k i/a,  where  p is  the  resistivity  of  the  plasma  and  O its 
P P 

conductivity.  From  Fig.  1.2  it  is  clear  chat  the  measured  voltage  SV 
is  related  to  these  other  quantities  by  the  equation 


inner 


6V,  + AV„  ^ „ - AV  , „ 

inner  W.F.,2  sheath, 2 


(k  i)/o  + AV„  „ , - AV  . . , 

p W.F.,1  sheath,! 


or,  upon  rearrangement, 
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where 
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offset 
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Now,  differentiating  Eqn.  (1.2)  (or  Eqn.  (1.3))  and  rearrangement  leads 
to  the  result 
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Figure  1.2.  Schematic  of  the  potential  distribution  along 
a line  joining  the  four  electrodes  of  a four- 
pin  probe. 
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in  which  it  is  understood  that  for  simplicity  “^^j^nner  denoted  by 
6V.  Since  k is  a constant,  the  conductivity  is  directly  proportional 
to  slope  of  the  current-voltage  characteristic  of  the  probe. 

The  probe  constant  must  be  known  to  at  least  the  accuracy  de- 

sired in  the  absolute  measurement  of  the  conductivity.  There  are  several 
approximate  ways  of  determining  k^.  Assuming  the  potential  distribution 
to  be  the  same  as  the  2-D  distribution  which  would  be  established  by  in- 
finitely long  probe  wires  of  vanishingly  small  diameter,  Brown  [6]  obtains 
the  following  formula  in  a straightforward  calculation: 

(!) 


Jln 


P.2-D 


ITS, 


(1.6) 


where  £ is  the  length  of  the  probe  wire  tips  and  a and  b are  geomet- 
rical parameters  defined  in  Fig.  1.1.  Formula  (1.6)  is  very  useful  for 
quick  order-of-magnitude  estimates  of  the  value  of  k^  for  any  given 
probe.  Somewhat  more  refined  estimates,  taking  into  account  the  finite 
diameter  of  the  electrodes,  could  be  obtained  by  solving  numerically,  or 
experimentally  with  resistance  paper,  the  two-dimensional  boundary -value 
problem  of  current  flow  between  two  parallel  cylindrical  electrodes  in  a 
uniform  medium  with  two  conducting  cylinders  in  between. 

In  principle  a similar  numerical  calculation  in  three  dimensions  is 
possible,  but  the  easiest  way  of  obtaining  k^  while  taking  account  of  the 
three-dimensionality  of  the  current  flow  around  the  pins  is  to  immerse  the 
probe  in  a liquid  electrolyte  (such  as  an  aqueous  KC£  solution).  If  the 
conductivity  of  this  solution  is  accurately  known  and  the  slope  of  the 
current-voltage  characteristic  can  be  measured  while  the  probe  is  in  the 
electrolyte,  then  k^  is  easily  calculated  using  Eqn.  (1.5).  This  was 
the  way  k^  was  in  fact  determined  for  the  probes  used  in  the  present  work. 
However,  several  difficulties  were  encountered  which  made  necessary  a fairly 
careful  and  elaborate  procedure  in  order  to  obtain  accurate  values  of  k^; 
the  details  of  this  probe  constant  determination  work  are  presented  in  Ap- 
pendix A.  It  is  worth  emphasizing  that  the  probe  constant  obtained  by  this 
method  is  valid  for  non-flowing,  stationary  conducting  media. 

In  practice,  then,  once  the  probe  constant  kp  has  been  determined, 
the  conductivity  of  the  test  plasma  may  be  found  by  Immersing  a four-pin 
probe  into  the  plasma,  exciting  the  probe  (l.e.,  passing  current  i at 
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some  suitable  frequency)  with  a power  source,  simultaneously  measuring  1 
and  the  inner-pin  voltage  difference  6V,  suitably  processing  the  i and 
6V  signals  so  as  to  obtain  the  slope  t and  calculating  a from 

the  relation  a = . The  time-response  of  the  probe  is  limited 

mainly  by  the  speed  with  which  the  electronic  signal  processing  and  record- 
ing can  be  carried  out.  The  spatial  resolution  depends,  of  course,  on  the 
size  of  the  actual  probe,  which  can  be  made  reasonably  small  (typical  di- 
mensions might  be  0.5  mm  diameter  wire  electrodes  of  2 or  3 mm  length  sepa- 
rated by,  say  1 mm).  There  is  a lower  limit  to  the  possible  probe  dimen- 
sions, however,  which  is  due  to  the  fact  that  the  plasma  sheaths  around  the 
inner  and  outer  pins  must  not  overlap  or  else  erroneous  measurements  of 

will  result.  In  addition,  it  should  be  noted  that  the  electrodes 
may  operate  at  sufficiently  high  temperatures  that  they  are  thermionically 
emitting  without  affecting  adversely  the  probe  reading.  Thermionic  emis- 
sion, in  fact,  simplifies  the  operation  of  the  probe  by  greatly  reducing 
the  voltage  drops  at  the  electrode-plasma  interface  and  permitting  large 
currents  to  be  passed  through  the  plasma  with  comparatively  small  applied 
voltages.  Operation  at  temperatures  of  thermionic  emission  also  elimi- 
nates the  need  for  elaborate  probe-cooling  systems,  thereby  simplifying 
the  construction  of  the  probes. 

1.4  The  Question  of  Accuracy  and  Flow  Effects 

In  the  introductory  discussion  of  the  four-pin  probe  in  the  previous 
section,  the  major  assumptions  have  been  that  (a)  the  probe's  i-6V  char- 
acteristic is  linear,  and  (b)  the  probe  constant  k^,  found  empirically 
in  non-flowing  electrolyte  solutions,  remains  valid  when  the  probe  is 
placed  in  a flowing  plasma.  These  assumptions  must  be  experimentally 
studied  if  a reliable  assessment  of  the  accuracy  and  suitability  of  this 
probe  method  under  MHD  conditions  is  to  be  made.  The  range  of  currents 
for  which  the  desired  linearity  holds  should  be  determined,  and  care  taken 
to  operate  the  probe  within  this  range.  The  possible  influence  of  "flow 
effects"  arises  when  one  considers  that  k^  is  a function  of  the  complete 
current  distribution  pattern  around  the  electrodes.  When  the  probe  is  in 
a uniform,  stationary  electrolyte,  a particular  current  pattern,  and  hence 
a certain  value  for  k^,  will  be  obtained.  But  when  the  probe  is  in  a 
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flowing  plasma,  there  are  several  effects  tending  to  produce  a non-uniform 
conductivity  field  around  the  pins,  such  as  thermal  boundary  layers  around 
the  pins  and  ceramic  support  which  might  lower  a in  regions  near  the 
cooler  surfaces,  or  convection  effects  tending  (in  loose  terminology)  to 
"blow  the  discharge  downstream".  Such  "flow  effects"  might  well  lead  to 
significantly  altered  (with  respect  to  the  stationary  medium  case)  current 
patterns  around  the  pins,  and  thus  to  an  effective  altered  value  of  k , 
or  in  other  words  to  an  indicated  conductivity  which  could  differ  signifi- 
cantly from  the  true  local  conductivity  of  the  plasma.  Just  how  much  the 
indicated  conductivity  depends  on  such  flow  effects  is  open  to  experimen- 
tal investigation.  For  convenience,  the  flow  effects  have  been  classified 
into  possible  "velocity  effects",  in  which  '^j^ndicated  ~ ^p  “ight 

change  as  the  magnitude  of  the  flow  velocity  changes  for  some  particular 
fixed  geometrical  orientation  of  the  pins  with  respect  to  the  velocity 
vector,  and  possible  "orientation  effects"  in  which  ^indicated  “^^ht 
change  as  the  angular  orientation  of  the  probe  relative  to  the  flow  changes 
while  the  magnitude  of  the  velocity  remains  constant. 


1 . 5 Specific  Goals  of  the  Research 

In  view  of  the  previous  discussion,  it  is  possible  now  to  state  briefly 
the  specific  goals  of  the  present  work.  The  first  major  goal  was  to  develop 
a suitable  electronics  probe  excitation  and  signal-processing  package  to 
measure  accurately  a.  ’ 


'indie  ^p  d(6V) 


on  both  a time-resolved  and  time- 


averaged  basis.  The  other  main  goals  are  to  provide  at  least  partial, 
hopefully  quantitative  answers  to  the  following  questions. 

• What  is  the  range  of  operating  parameters  (if  any)  for  which  the  i-6V 
characteristics  of  the  four-pin  probe  is  linear? 

• For  a fixed  orientation  of  the  four-pin  probe,  and  fixed  a,  does 
Olndic  depend,  and  if  so  by  how  much,  on  the  magnitude  of  the  flow 
velocity? 

• For  a fixed  flow  velocity,  does  *^jjjdlc  the  angular  orien- 

tation of  the  pins  is  changed  relative  to  the  flow,  and  if  so,  by  how 
much? 

• To  what  accuracy  or  degree  of  confidence  can  one  assert  that  '^iudic 
as  measured  by  a four-pin  probe  is  equal  to  the  true  local  electrical 
conductivity  in  flowing  MHD  plasmas? 
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1.6  Summary  of  What  Follows,  and  of  Conclusions 


Chapter  2 of  this  report  describes  the  experimental  apparatus  other 
than  the  electronic  instrumentation:  the  arcjet  facility  and  flow  train 
for  generating  the  seeded  argon  plasma,  as  well  as  construction  details 
of  the  four-pin  probes. 

The  experiments  to  investigate  the  i-6V  characteristics  of  four-pin 
probes,  and  establish  the  range  of  linearity,  are  discussed  in  Chapter  3. 

The  results  of  these  experiments  were  important  in  motivating  the  final 
electronics  signal-processing  system,  which  is  the  main  topic  of  Chapter  4. 

Once  the  final  electronics  configuration  was  implemented,  a number 
of  experiments  were  done  using  this  system  to  make  both  time-averaged  and 
time-resolved  measurements  of  O.  The  time-averaged  measurements  were 
mainly  to  test  the  sensitivity  of  the  Indicated  conductivity  to  flow  ef- 
fects (velocity  and  orientation);  these  results  are  presented  in  Chapter  5. 
Chapter  6 presents  some  sample  time-resolved  data  and  the  statistical  in- 
formation about  the  conductivity  fluctuations  which  was  deduced  from  such 
data. 

The  main  conclusions  which  we  have  drawn  about  four-pin  probes  in 
the  course  of  this  research  are  as  follows: 

• The  i-6V  characteristics  of  the  four-pin  probes  used  remained  linear 
if  the  voltage  applied  to  the  outer  pins  was  maintained  below  approximately 
1 volt.  Depending  on  the  particular  probe  geometry,  the  corresponding  cur- 
rent densities  at  the  outer  pin  surfaces  were  less  than  approximately 

2 2 

0.1  arap/cm  to  0.2  amp/cm  . The  criterion  V ^ <1  volt  seemed  to 

^ outer 

be  the  most  practical  and  reliable  rule  of  thumb  for  ensuring  linearity. 

• Even  when  no  current  passes  between  the  outer  electrodes,  a finite 
and  fluctuating  voltage  difference  ‘^'^offset  appears  between  the  inner 
pins;  this  fact  requires  that  the  probe  not  be  operated  with  d.c.  current 
if  erroneous  results  are  to  be  avoided. 

• When  excited  at  frequencies  well  above  the  highest  naturally-occurring 
conductivity  fluctuation  frequencies  in  the  flow,  the  four-pin  probe  system 
described  in  this  report  appears  to  yield  reliable  and  accurate  time-resolved 
data  on  conductivity  fluctuations  (accurate  in  the  sense  that  all  statisti- 
cal parameters  characterizing  the  fluctuations  are  accurate,  with  the  excep- 
tion of  the  mean  value  of  o) . 
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I • With  the  method  of  immersion  in  electrolyte  solutions  to  calibrate 

the  probes,  can  be  found  only  to  a precision  of  approximately  5%  to 

10%  (depending  on  the  probe  used),  which  limits  the  probe's  absolute  accu- 
racy even  in  the  absence  of  flow  effects. 

• For  a fixed  probe  orientation,  the  indicated  conductivity 

i decreases  as  the  magnitude  of  the  flow  velocity  of  the  plasma  increases. 

t 

: The  curve  of  “ presented  in  Fig.  5.3  shows  a drop  of  ap- 

j proximately  50%  in  ^^jjdic  increases  from  about  10  m/sec  to  100  m/sec. 

I • For  a fixed  flow  velocity,  *^indic  appears  to  decrease  as  the  flow 

direction  changes  from  "end-on"  towards  a sidewise  orientation  relative  to 
^ the  pins  (cf.  Fig.  5.6). 

' • The  present  work  has  not  studied  possible  effects  of  the  size  and 

geometry  of  the  probe  on  has  restricted  attention  mainly  to 

a particular  probe  in  the  "end-on"  configuration;  nevertheless,  it  has 
shown  that  with  careful,  detailed  measurements  it  is  possible,  for  any  par- 
r ticular  probe  and  orientation,  to  produce  an  empirical  curve  of  ®indic^^u=o 

vs.  u which  can  be  used  as  a flow-effect  calibration  curve,  i.e.,  to 
provide  a correction  to  *^j^ndic  '^hlch  would  yield  an  accurate  measure  of 
the  true  local  conductivity  at  any  velocity. 

• In  the  absence  of  such  detailed  flow-effect  calibration  measurements 

as  are  mentioned  in  the  preceding  conclusion,  it  would  be  unwarranted  (based 
on  the  data  obtained  in  this  work)  to  claim  high  absolute  accuracy  for  any 
particular  four-pin  probe  used  in  flowing  MHD  plasmas,  since  the  indicated 
conductivity  may  be  as  much  as  a factor  of  two  or  three  lower  than  the 
actual  conductivity. 

• Care  should  be  exercised  in  Interpreting  conductivity  profile  data  ob- 

j tained  by  four-pin  probes  (such  as  some  of  the  results  of  Brown  [6])  in 

1 boundary  layers  where  both  a and  the  velocity  are  varying. 


Chapter  2 
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DESCRIPTION  OF  THE  EXPERIMENTAL  APPARATUS 

2 . 1 Probe  Construction  Details 

All  of  the  four-pin  probes  used  in  this  work  were  constructed  with 
tungsten  wire  electrodes  supported  by  double-bore  alumina  tubing.  No 
artificial  cooling  was  used,  since  these  materials  can  survive  the  NaK- 
seeded  argon  plasma  at  around  2000°K  without  melting  or  oxydizing. 

There  were  two  basic  variations  in  the  actual  construction  of  the 
probes.  In  the  simplest  version,  which  was  employed  as  a fixed-position 
probe  to  monitor  conductivity  changes  in  the  plenum,  the  two  ceramic  tubes 
were  inserted  through  holes  drilled  into  a brass  pipe  plug  which  was  screwed 
into  the  water-cooled  plenum  wall,  as  shown  in  Fig.  2.1.  The  tungsten 
wires  were  typically  0.020  or  0.030  inch  diameter,  and  the  exposed  tips 
were  located  close  to  the  centerline  in  the  low-velocity  region  in  the 
plenum. 

The  other  type  of  probe  construction  was  used  for  mobile  probes  which 
were  mounted  on  a mechanical  traversing  mechanism  situated  outside  the 
plenum  of  the  plasma  generator  facility.  The  ceramic  tubes  for  these  mobile 
probes  were  inserted  into  5/8"  diameter  graphite  rods  appropriately  drilled 
and  slotted,  as  Indicated  in  Fig.  2.2.  The  alumina  tubes  were  potted  into 
the  holes  and  slots  in  the  graphite  with  Ceramabond  503  adhesive.  In  the 
case  of  long  probes  employing  '.hln  ceramic  tubing,  it  was  necessary  to 
wrap  the  tubing  with  thin  molybdenum  wire,  as  shown  in  Fig.  2.2,  to  keep 
the  two  tubes  from  separating  due  to  thermal  warping;  some  probes  were  also 
reinforced  by  securing  1/16  inch  diameter  tungsten  rods  along  the  outside 
of  the  alumina  tubes  to  reduce  lateral  warping.  While  liberal  use  of 
Ceramabond  adhesive  was  made  to  hold  the  tubes  together  and  coat  the  ex- 
posed metal  surfaces  of  the  moly  wire  and  tungsten  reinforcing  rods,  care 
was  taken  not  to  place  any  Ceramabond  in  the  region  of  the  exposed  tungsten 
electrode  pins.  The  exposed  pins  were  cleaned  with  a dental  sandblaster 
before  use,  and  appeared  clean  and  bright  after  each  test. 
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Water-cooled 
Plenum  Wall 


Figure  2.1.  Construction  of  fixed  probes.  Flow  direction  is  normal 
to  the  plane  of  the  page. 


Plasma  Flow 


5/8-inch  diameter 
graphite  rode 


Figure  2.2.  Construction  details  of  mobile  probes.  The  alumina 
tubing  was  potted  into  the  graphite  with  Ceramabond 
503  high-temperature  cement,  which  is  not  shown  in 
this  sketch. 
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2 , 2 Plasma-Production  Apparatus  and  Flow  Path 


Argon  seeded  with  NaK  was  employed  as  the  working  gas  in  all  experi- 
ments. The  argon  was  heated  by  passing  the  gas  through  the  swirl-stabilized 
discharge  in  a Thermal  Dynamics,  Inc.  Model  F-5000  arc.  heater.  The  argon 
flowrate  was  normally  set  at  13.6  gm/sec,  and  the  are  operated  at  about 
700  amperes  d.c.  with  a voltage  drop  of  about  90  volts;  the  temperature  of 
the  gas  usually  stayed  in  the  range  ITOO^K  to  2300®K,  depending  on  the  arc 
current  setting.  To  smooth  out  the  arcjet  operation,  a large,  water-cooled, 
iron-core  conventional  magnet  (normally  used  in  combustion  MHD  research  in 
our  laboratory)  was  used  in  series  with  the  arcjet.  Its  large  inductance 
was  very  effective  in  filtering  out  the  360  Hz  ripple  in  the  arcjet  cur- 
rent (which  led  to  360  Hz  conductivity  fluctuations  observed  in  early  tests) 
arising  from  the  three-phase  full-wave  rectified  power  source  employed  in 
this  facility.  Further  details  of  this  aspect  of  the  experiments  are  given 
in  Appendix  B. 

After  leaving  the  arc  itself,  the  heated  argon  flows  into  a cylindri- 
cal water-cooled  plenum  chamber,  as  indicated  in  Fig.  2.3.  The  molybdenum 
liner  shown  in  this  figure  acts  as  a radiation  shield  to  decrease  the  heat 
transfer  rate  to  the  water-cooled  plenum  wall.  At  the  downstream  end  of 
the  plenum  the  gas  flows  through  a slowly  converging  nozzle  made  of  boron 
nitride,  and  exhausts  into  the  atmosphere. 

Seeding  is  accomplished  by  injecting  the  liquid  metal  NaK  (sodium- 
potassium  mixture,  78%  potassium  by  weight)  through  a molybdenum  tube  into 
the  plenum.  As  the  NaK  enters  the  heated  tube  it  vaporizes  and  emerges 
from  a single  0.020  inch  diameter  hole  facing  downstream  (it  was  found  that 
when  this  hole  faced  upstream,  some  seed  presumably  entered  the  arc  region 
and  interfered  with  the  stable  operation  of  the  arcjet).  Further  details 
on  this  seed-injection  system  may  be  found  in  Reference  [7].  To  avoid  the 
condensation  of  liquid  sodium  and  potassium  metal  on  the  boron  nitride  noz- 
zle walls,  a stainless  steel  insert  was  placed  between  the  nozzle  and  the 
cooled  copper  end  flange.  This  insert  was  designed  with  large  thermal 
resistance  and  allowed  the  nozzle  to  run  at  temperatures  above  1000°K, 
thus  avoiding  NaK  condensation. 


2.3  Mobile  Probe  Mount  and  Technique  for  Varying  Velocity 


Since  one  of  the  main  goals  of  the  present  work  was  to  investigate 
the  effects  of  flow  velocity  and  orientation  on  the  indicated  conductivity 
reading,  a method  was  needed  for  easily  varying  both  the  flow  velocity  over 
the  probe  and  the  probe  orientation  with  respect  to  the  flow  direction, 
while  leaving  the  local  electrical  conductivity  of  the  plasma  approximately 
constant.  It  was  impractical  to  attempt  varying  the  flow  velocity  by 
varying  the  argon  flowrate  through  the  arcjet  since,  in  order  to  keep  the 
conductivity  constant,  both  the  arcjet  current  (and  hence  the  gas  tempera- 
ture) and  the  seeding  rate  would  have  to  be  carefully  readjusted.  It  was 
therefore  decided  to  employ  a slowly  converging  nozzle  at  the  plenum  exit; 
the  flow  velocity  at  each  point  along  the  axis  of  the  nozzle  would  then  be 
different  (with  the  conductivity  remaining  essentially  constant,  since  the 
arcjet  current  and  seed  rate  were  in  steady  state).  A mechanical  mount 
was  constructed  using  a translating  rotary  milling  head  to  insert  and  with- 
draw a mobile  four-pin  probe  along  the  nozzle  axis  in  the  region  of  varying 
velocity,  as  shown  in  Fig.  2.4.  'Ihe  axial  position  of  this  probe  was  moni- 
tored continuously  with  a linear  position  transducer.  The  calibration  pro- 
cedure by  which  this  axial  position  was  related  to  the  corresponding  flow 
velocity  is  discussed  in  Appendix  C.  Using  this  system,  the  flow  velocity 
varied  from  about  100  m/sec,  with  the  mobile  probe  fully  withdrawn  out  to 
nearly  the  exit  plane  of  the  nozzle,  down  to  about  10  m/sec,  with  the  probe 
fully  Inserted  into  the  plenum  region. 

By  placing  the  axis  of  rotation  of  the  rotary  milling  head  directly 
beneath  the  exposed  tips  of  the  mobile  probe  and  translating  the  probe 
tips  to  a point  just  inside  the  nozzle  exit  plane,  the  angular  orientation 
of  the  pins  could  be  continuously  varied  by  rotating  the  milling  head. 

The  probe  tips  remain  at  the  same  radial  and  axial  location  near  the  noz- 
zoe  exit  during  this  process  so  both  the  conductivity  and  the  magnitude 
of  the  flow  are  unchanged.  The  range  of  variation  of  the  orientation  angle 
6 (see  Fig.  2.4b)  was  about ±60°;  a full  ±90°  variation  was  not  possible, 
since  this  would  have  required  that  the  probe  tips  lie  slightly  outside  the 
exit  plane  of  the  nozzle  and  it  was  found  in  early  tests  that  in  this  case 
when  6 = 90°  (l.e.,  when  the  probe  was  Inserted  from  the  side  normal  to 
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the  flow) , some  atmospheric  oxygen  was  entrained  along  the  downstream 
side  of  the  ceramic  support  tubes  and  reached  the  exposed  tungsten  pins, 
which  quickly  oxidized. 

An  additional  degree  of  freedom  in  the  translating  mobile  probe  mount 
permitted  traversing  the  probe  radially  across  the  stream  and  into  the 
boundary  layers.  This  was  useful  in  measuring  the  severity  of  the  radial 
gradients  (or  in  measuring  the  radial  profiles)  of  the  conductivity. 
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Chapter  3 

INVESTIGATION  OF  i-6V  CHARACTERISTICS  OF  FOUR-PIN  PROBES 


3.1  Preliminary  Instrumentation  for  Studying  Probe  Characteristics 

The  first  system  for  varying  the  current  passing  between  the  outer 
electrodes  of  the  probe,  and  for  generating  directly  plots  of  this  current 
(i)  vs.  the  inner  pin  potential  difference  (6V),  is  the  simple  battery- 
driven  circuit  and  x-y  plotter  shown  schematically  in  Fig.  3.1.  With 
the  probe  immersed  in  the  seeded  argon  plasma,  the  potentiometer  was  slowly 
turned  by  hand  from  one  extreme  to  the  other  over  a time  period  of  the 
order  of  one  or  two  minutes,  and  the  signals  6V,  the  inner  pin  voltage 
difference,  and  V , ^ = iR  . ^ (used  as  a measure  of  the  current)  were 

applied  directly  to  the  x-  and  y-axes,  respectively  of  the  x-y  plotter, 
whose  differential  input  terminals  were  floating  off  ground.  (In  all  ex- 
periments with  four-pin  probes,  care  was  taken  not  to  ground  any  point  in 
the  probe-excitation  circuitry,  and  all  instruments  to  which  the  ^g^unt 
and  6V  signals  were  applied  had  to  be  true  floating,  differential-input 
devices.  These  precautions  were  necessary,  because  the  probe  pins  them- 
selves were  immersed  in  a plasma  which  was  in  good  electrical  contact  with 
the  grounded  anode  of  the  arcjet  facility;  thus  any  other  grounded  point 
in  the  probe  power  supply  would  provide  a possible  alternative  path  for 
the  probe  current,  part  or  all  of  which  might  bypass  the  shunt  resistor, 
causing  erroneous  measurements.) 

A reproduction  of  a typical  1-6V  plot  obtained  in  this  way  is  shown 
in  Fig.  3.2.  While  the  general  appearance  of  this  1-6V  characteristic 
indicates  that  the  probe  is  operating  approximately  as  expected,  there  are 
several  deviations  from  ideal  behavior.  First,  the  trace  is  very  Jagged, 
especially  near  the  extremes.  This  behavior  could  be  due  in  part  to  uneven- 
ness in  the  manual  rotation  of  the  potentiometer,  and  in  part  due  to  either 
fluctuations  in  plasma  properties  during  the  long  sweep  time  or  to  some 
sort  of  arcing  at  the  probe  pins.  Secondly,  the  overall  shape  of  the  i-6V 
trace  is  slightly  "S-shaped"  and  nonlinear,  although  the  portion  near  the 
origin  appears  linear  over  some  range  of  currents.  The  slope  of  the  trace 


22 


Figure  3.1.  Initial  circuit  for  generating  i-6V  characteristics 
of  4-pin  probes. 
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Probe  current 


Reproduction  of  typical  i_(Sv  characteristic  obtained  with  the  circuit  of 
Figure  3.1. 
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taken  near  the  origin  does  lead  to  a value  of  conductivity  consistent  with 
the  plasma  temperature  (as  measured  only  approximately  by  thermocouples) 
and  seed  fraction  present  in  the  experiment.  Thirdly,  as  the  current  is 
raised  and  then  lowered,  the  plotter  pen  does  not  retrace  the  original 
curve;  indeed  there  even  appears  to  be  a sort  of  "hysteresis  effect"  as 
the  current  is  taken  through  a complete  cycle.  Additionally,  the  time 
(1-2  minutes)  required  to  manually  sweep  the  current  is  so  long  that  appre- 
ciable changes  in  the  conductivity  may  well  occur  as  the  i-6v  plot  is 
being  made. 

3 . 2 Revised  Measurement  Circuitry  and  Instrumentation  Setup 

The  difficulties  mentioned  in  the  previous  paragraph  motivated  several 
revisions  in  the  probe-excitation  circuitry  and  recording  instrumentation. 
To  replace  the  battery  excitation,  a completely  floating  power  supply  was 
constructed,  based  on  a power  operational  amplifier.  Any  desired  excita- 
tion current  waveform,  from  pure  d.c.  to  a.c.  at  frequencies  up  to  near 
100  kHz,  could  be  obtained  by  applying  a corresponding  low-power,  floating 
excitation  voltage  to  the  input  of  the  power  operational  amplifier.  The 
complete  schematic  and  further  description  of  this  floating  power  amplifier 
appears  in  Appendix  D. 

For  purposes  of  studying  the  time-averaged  i-6V  characteristics,  an 
excitation  source  was  needed  which  was  faster  and  smoother  than  the  manual 
sweeping  of  a potentiometer  (as  described  in  Section  3.1)  while  still  slow 
enough  so  that  the  x-y  plotter  pen  could  accurately  track  the  signals. 

For  this  purpose,  a small  battery-powered  function  generator  was  construc- 
ted from  a kit  (see  Appendix  D for  further  details);  this  device  produces 
a sinusoidal  or  triangular  waveform  at  frequencies  from  0.1  Hz  to  100  kHz, 
and  was  operated  at  0.1  Hz  for  purposes  of  obtaining  1-6V  characteristics 
on  the  X-Y  plotter. 

In  order  to  "clean  up"  the  noisy,  fluctuating  traces  typified  by  Fig. 
3.2,  the  6V  and  signals  are  passed  through  low-pass  filters 

before  being  sent  to  the  plotter.  With  the  excitation  frequency  set  at 
0.1  Hz,  the  cutoff  frequency  of  these  commercial  filters  (Rockland  Systems 
dbrp.  Model  452  Dual  Hl/Lo  Filters)  was  typically  set  at  0.2  Hz.  Since  the 
input  terminals  on  these  filters  are  single-ended  (i.e.,  one  input  terminal 
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is  grounded) , and  It  was  desired  to  maintain  the  entire  measuring  circuit 
floating,  it  was  necessary  to  pass  the  6V  and  signals  through 

unity-gain  differential  amplifiers,  which  act  as  isolation  stages,  before 
Inputting  these  signals  into  the  filters.  With  all  of  these  changes,  the 
final  experimental  electronics  configuration  for  obtaining  "high  quality" 
1-6V  plots  from  the  four-pin  probe  is  as  shown  in  Fig.  3.3. 

3.3  Results  of  Studies  of  i-6V  Characteristics 

Using  the  experimental  setup  of  Fig.  3.3,  i-6V  characteristics  of 
two  different  probes  were  recorded  after  varying  (and  recording)  the  ampli- 
tude of  the  0.1  Hz  voltage  applied  to  the  outer  pins  of  the  probe, 

Vouter  max’  varying  the  plasma  temperature  (and  thus  the  magnitude  of 

the  conductivity).  In  several  of  the  runs,  the  seeded  argon  was  diluted 
with  3%  by  volume  of  nitrogen  to  reduce  the  possibility  of  electron  tem- 
perature nonequilibriura  in  regions  of  locally  high  current  density.  For 
the  various  run  conditions,  simultaneous  plots  of  i vs.  V and  of 

i vs.  6V  were  recorded.  The  main  results  of  these  experiments  will 
now  be  discussed. 

It  was  found  that  the  i-6V  characteristics  were  always  linear  so 

long  as  V ^ was  kept  below  about  1.0  volts.  Also,  for  values  of 

V ^ less  than  roughly  0.5  volts,  the  quality  of  the  i-6V  plots 

outer,  max  ^ y j f 

seemed  to  deteriorate,  presumably  due  to  decreased  signal-to-noise  ratios. 

Sample  1-6V  plots  for  varying  V are  reproduced  in  Figs.  3.4a, 

^ ° outer, max  ^ ® 

3.4b,  and  3.4c,  in  which  V ^ had  values  of  6 volts,  0.17  volts, 

outer, max 

and  0.71  volts,  respectively.  The  best  quality,  most  linear  i-6V  char- 
acteristics were  obtained,  no  matter  which  probe  was  used  and  no  matter 

what  the  numerical  value  of  o,  for  V =0.75  volts,  and  Fig. 

outer, max  ° 

3.4c  is  a typical  sample  characteristic.  It  is  clear  that  the  low-pass 
filters  smooth  out  the  small  fluctuations  such  as  those  in  Fig.  3.2,  and 
the  "hysteresis  effect"  is  minimized  or  absent  entirely.  The  values  of 
a obtained  by  measuring  the  slopes  of  the  characteristics  and  multiplying 
by  the  probe  constant  obtained  using  the  formula  (1.6)  always  lie  within 
the  (rather  wide)  range  of  expected  values  of  o based  on  the  relatively 
inaccurate  thermocouple-determined  temperature.  The  behavior  of  the  i-6V 
characteristics  was  unchanged,  so  far  as  we  were  able  to  judge,  by  the 
presence  of  nitrogen  in  the  plasma. 
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Figure  3.4.  Sample  i-6V  characteristics  of  four-pin  probes  for 
several  values  of  from  test  of  10/29/76. 
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It  should  be  noted  that  the  experimental  plots  of  i vs.  6V 

outer 

were  not,  except  in  a few  cases,  symmetrical  with  respect  to  the  origin, 

nor  were  they  always  linear,  even  when  Iv  1-1  volt.  The  inter- 

' outer, max' 

pretation  of  these  ^"^'^outer  felt  to  be  somewhat  Irrelevant 

to  the  main  thrust  of  the  present  work,  however,  especially  since  even  in 

cases  where  the  ^“'^^outer  linear,  the  corresponding  "more 

important"  i-6V.  curve  was  linear,  provided  |v  1-1  volt, 

inner  ^ ' outer, max' 

Before  carrying  out  these  experiments,  we  had  expected  to  be  able  to 

formulate  a rule-of-thumb  for  obtaining  linear  i-6V  characteristics  in 

terms  of  the  current  density  at  the  surface  of  the  outer  probe  electrodes. 

However,  the  results  showed  that  among  all  "good  quality",  linear  i-6V 

characteristics  which  were  recorded,  this  current  density  varied  over  the 
2 2 

range  0.05  amps/cm  to  0.25  amps/cm  , depending  on  the  probe  size  and 

on  the  plasma  temperature.  The  most  useful  guide  to  ensuring  a linear 

characteristic  appeared  to  be,  as  has  been  stated  already,  to  maintain 

V ^ in  the  range  0.5  volts  to  1.0  volts,  with  the  cleanest  data 

outer, max  ° ’ 

found  for  V ^ =0.75  volts;  then,  so  long  as  the  characteristic  was 

outer, max  ’ ® 

linear,  the  current  density  increased  proportionately  with  the  conductivity. 

Perhaps  the  most  important  fact  to  emerge  from  studying  the  experimen- 
tal i-6V  plots,  at  least  insofar  as  the  subsequent  course  of  the  experi- 
ments was  concerned,  was  the  fact  that  even  when  the  probe  current  i goes 
to  zero,  there  is  still  a fairly  large  voltage  difference  between  the  inner 
pins,  denoted  by  '^^gffset  found  graphically  by  locating 

where  the  i-6V  trace  intercepts  the  horizontal  (6V)  axis.  Experimen- 
tally, for  a given  probe,  the  value  of  6V  tended  to  fluctuate 

offset 

slowly  around  some  mean  value  (for  example,  for  the  probe  used  to  obtain 
Fig.  3.4c,  "^^offset  ~ mV).  The  particular  mean  value, 

6Vof fset>  however,  seemed  to  vary  randomly  from  one  probe  to  the  next, 
usually  falling  within  the  range  -50  mV  < fggj-  ~ '*■5®  precise 

origin  of  this  6V  is  not  known.  In  any  case,  its  magnitude  is  at 

most  a few  percent  of  the  floating  potential  (p^  = (kT/e) £.n(m^/m^) 

(cf.  Reference  [8],  p.  132,  Eqn.  (3.5)),  which  for  the  experimental  condi- 
tions of  this  work  has  a value  of  approximately  1.0  volts.  As  has  been 
suggested  in  Section  1.3  (cf.  Fig.  1.2  and  Eqn.  (1.4)),  some  plausible 

mechanisms  giving  rise  to  a nonzero  6V  are  differences  in  sheath 

offset 
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voltage  drops  (due  perhsps  to  slight  differences  in  probe  tip  temperatures) 
and/or  differences  in  the  work  function  voltage  drops.  The  effective  work 
function  at  the  surface  depends  on  the  contact  potential  between  the  tung- 
sten electrode  and  the  partial  coating  of  potassium  seed  atoms.  If  the 
seed  atom  coverage  of  one  of  the  two  inner  floating  electrodes  differs 
slightly  from  the  coverage  on  the  other  pin,  then  there  could  be  differ- 
ences in  the  contact  potentials  which  would  show  up  as  ‘^'^offset 
potential  difference  between  the  two  inner  pins  is  measured. 

Whatever  the  mechanism  leading  to  6V  , its  existence  had  to 

® offset 

be  taken  into  account  in  designing  subsequent  signal-processing  electronics. 

What  is  desired  to  obtain  0 is  the  slope  ■ of  the  i-6v  character- 

u Co  V/ 

Istic.  If  the  characteristic  is  linear  and  passes  through  the  origin  of 
coordinates  (refer  to  Fig.  3.5),  then  this  slope  would  be  the  quotient  of 
the  amplitudes  of  the  two  a.c.  signals,  which  is  easy  to  obtain  continu- 
ously using  an  electronic  analog  divider.  However,  if  (as  is  in  fact  the 
case)  the  characteristic  is  linear  but  does  not  pass  through  the  origin, 
then  the  simple  quotient  of  the  maximum  i and  6V  signals  gives  a com- 
pletely erroneous  value  for  the  slope  (see  Fig.  3.5),  and  hence  for  the 
conductivity.  The  implication  from  all  this  is  that  a four-pin  probe 
should  not  be  operated  in  the  d.c.  mode  (i.e.,  with  d.c.  current  passing 
between  the  outer  pins),  but  rather  in  the  a.c.  mode,  in  which  the  i-6V 
characteristic  is  being  continuously  swept  out.  Then,  if  the  '^'^offset 
voltage  can  be  removed  electronically,  the  quotient  of  the  amplitudes  of 
the  1 and  6V  signals  is  equal  to  the  slope  of  the  characteristic,  and 
hence  is  proportional  to  the  conductivity.  The  details  of  how  this  was 
accomplished  in  the  present  work  are  presented  in  the  following  chapter. 
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(a)  6 Voffs^t  = 0 0*^  removed  electronically  by  signal 

processing 
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Figure  3.5.  Illustrating  possible  error  in  obtaining  true 
slope  of  I-6V  characteristics  as  the  quotient 
Imax/'5Vmax  presence  of  non-zero 

'^'/offset’ 
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Chapter  4 


FINAL  SIGNAL-PROCESSING  SYSTEM 

4 . 1 Use  of  High-Pass  Filters  and  Audio-Frequency  Excitation 

In  order  to  obtain  an  1-6V  characteristic  which  was  linear  and 
always  passed  through  the  origin,  it  was  necessary  to  filter  electroni- 
cally the  i and  6V  signals.  Since  the  offset  signals  were  observed 
to  vary  at  frequencies  less  than  one  or  two  kilohertz,  it  was  decided  to 
power  the  outer  pins  of  the  probe  at  a much  higher  frequency  (well  into 
the  audio-frequency  range)  and  to  employ  a high-pass  filter  adjusted  to 
attenuate  the  offset  signal,  but  to  pass  the  excitation  signal.  A typical 
value  used  for  the  excitation  frequency  was  20  kHz,  in  which  case  the  high- 
pass  filter's  cutoff  frequency  would  be  set  around  15  kHz. 

The  instrumentation  setup,  shown  below  in  Fig.  4.1,  is  similar  in 
appearance  to  the  low-frequency  setup  of  Fig.  3.3  previously  described, 
except  that  now  the  source  frequency  was  raised  from  0.1  Hz  to  around 
20  kHz,  and  the  low-pass  filters  with  cutoff  frequencies  around  0.2  Hz 
were  switched  to  high-pass  filters  with  cutoffs  of  10  kHz  to  15  kHz. 
(Actually  the  same  instruments,  Rockland  Model  452  Dual  Hi/Lo  Filters, 
which  can  be  operated  in  either  the  high-pass  or  low-pass  mode  and  have 
adjustable  cutoff  frequencies,  were  used  in  both  setups.)  The  differential 
amplifiers  again  play  the  role  of  isolating  the  two  signals  from  ground, 
and  were  now  operated  with  sufficient  gain  to  raise  the  signal  levels  from 
a few  tens  of  millivolts  to  a few  volts  before  Inputting  them  into  the 
high-pass  filters.  Since  the  X-Y  plotter  clearly  cannot  plot  out  char- 
acteristics at  20  kHz,  an  oscilloscope  was  used  as  the  X-Y  plotting  de- 
vice, with  the  amplified  and  filtered  current  signal  applied  to  its  verti- 
cal deflection  plates  and  the  corresponding  6V  signal  applied  to  the 
horizontal  deflection  plates. 

A typical  photograph  of  the  resulting  i-6V  characteristics  is  re- 
produced in  Fig.  4.2.  Note  that  the  plots  are  all  very  linear  and  do  pass 
through  the  origin.  Also,  there  is  no  noticeable  phase  shift  between  i 
and  6V,  since  the  plots  are  linear  and  not  elliptical.  During  the 
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Figure  4.2.  Typical  i-6V  plot  obtained  with  the 

measuring  circuit  of  Figure  4.1.  Data 
from  run  133  of  2/11/77.  The  minimum 
and  maximum  slopes  correspond  to 
conductivities  of  10.0  mho/m  and 
26.8  mho/m,  respectively.  Exposure 
time  = 1/60  sec. 
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exposure  time  for  this  photograph  (1/60  sec),  many  complete  i-6V  plots 
were  swept  out,  and,  while  the  Individual  plots  are  linear,  their  slopes 
vary  over  some  range.  As  is  shown  in  more  mathematical  detail  in  Appen- 
dix E,  these  variations  in  slope  actually  represent  and  are  proportional 
to  the  fluctuations  in  the  electrical  conductivity  in  the  plasma,  even 
though  such  fluctuations  may  have  frequencies  well  below  the  cutoff  fre- 
quencies of  the  high-pass  filters.  The  basic  reason  for  this  at-first^ 
confusing  result  (viz,  how  low-frequency  conductivity  fluctuations  can 
"get  through"  the  high-pass  filters)  is  that  the  conductivity  fluctuations 
affect  only  the  amplitudes  of  the  high-frequency  "carrier  signals"  i and 
6V,  which  pass  unaltered  (in  either  frequency  or  amplitude)  through  the 
filters.  Thus  the  amplitudes  of  the  signals  emerging  from  the  filters  can 
and  do  still  vary  at  the  frequencies  of  the  (slower)  conductivity  varia- 
tions, which  consequently  are  seen  as  changes  in  the  slope  of  the  i-6V 
characteristics,  such  as  those  of  Fig.  4.2. 


4 . 2 Remaining  Signal-Processing  Components 

As  a result  of  the  signal  processing  described  in  the  previous  sec- 
tion, the  signals  emerging  from  the  two  high-pass  filters  have  the  form 


V , .(t)  = V , ^ (t)  sin  w t = R , I (t)  sin  u)  t 

shunt  shunt, max  s shunt  max  s 


(4.1) 


6V(t)  = 6V  (t)  sin  o)  t , 

max  s ’ 


where  is  the  excitation  frequency  20  kHz).  Since  the  plot  of 

Vshunt^^^  vs.  6V(t)  has  the  form  of  straight  lines  passing  through  the 

origin  of  coordinates  (see  Fig.  4.2),  the  slope  of  the  characteristic 

Vshunt  is  proportional  to  the  conductivity,  is  just  equal 

to  the  quotient  V , (t)/6V  (t),  so  we  may  write 

shunt, max  max 

^ . (t) 

0(t)  = (const.)  — gy — ^-T^) (^.2) 

max' 

Thus,  in  order  to  continuously  monitor  o(t),  one  needs  to  obtain  signals 
directly  proportional  to  the  amplitudes 

input  these  into  the  numerator  and  denominator  input  terminals  of  an 
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analog  divider  circuit,  whose  output  will  then  be  a signal  proportional 
to  the  time-varying  conductivity,  o(t). 

Signals  proportional  to  V , ^ (t)  and  6V  (t)  were  obtained 

® r r shunt, max  max 

by  first  passing  both  the  signals  represented  by  Eqn.  (4.1)  through  pre- 
cision linear  full-wave  rectifier  circuits  (see  Appendix  D for  circuit 
details),  and  then  sending  the  rectified  signals  into  low-pass  filters 
with  cutoff  frequencies  well  below  the  excitation  frequency  but  above 

the  highest  expected  frequency  of  the  conductivity  fluctuations.  The  out- 
puts of  these  latter  low-pass  filters  consisted  of  time-varying  positive- 
valued voltages  proportional  to  the  amplitudes 

dVmax^t).  These  latter  signals  were  subsequently  amplified  by  variable- 
gain  amplifiers  whose  gains  were  selected  for  each  experiment  to  give  opti- 
mum input  values  for  the  numerator  and  denominator  of  an  Analog  Devices 
Model  436A  analog  divider.  (This  particular  device  requires  that  at  all 
times  the  numerator  input  signal  be  less  than  or  equal  to  the  denominator 
input  signal,  and  that  the  denominator  signal  lie  within  the  range  0.1  volt 
to  10  volts.  The  variable-gain  amplifiers  allowed  the  operator  to  adjust 
the  signal  levels  to  comply  with  these  conditions  before  taking  data.) 

A schematic  diagram  showing  the  effects  of  the  rectifiers,  low-pass  fil- 
ters, and  variable-gain  amplifiers  on  a sample  waveform  emerging  from  the 
high-pass  filter  is  shown  in  Fig.  4.3  below. 

The  voltage  emerging  from  the  analog  divider  is  directly  proportional 
to  the  electrical  conductivity  of  the  plasma  as  a function  of  time,  and 
provides  a time-resolved  measurement  of  C'(t).  To  record  this  time-resolved 
conductivity  data,  the  signal  was  sent  directly  to  an  analog  tape  recorder, 
and  also  to  an  oscilloscope  for  visual  monitoring.  For  experiments  in 
v/hich  the  time— averaged  conductivity  was  the  variable  of  interest,  this 
same  signal  was  sent  to  a final  low- pass  filter  which  smoothed  out  the 
fluctuations,  and  thereafter  to  a stripchart  recorder  which  gave  a contin- 
uous visual  display  of  the  mean  conductivity.  The  procedure  for  calibrat- 
ing this  system,  taking  into  account  the  gain  factors  of  the  several 
electronic  signal-processing  components  and  the  probe  calibration  constant  , 

kp  (cf.  Section  1.3)  is  discussed  in  the  following  two  sections  of  this  ; 

chapter.  First,  however,  we  present  in  Fig.  4.4  a block  diagram  of  the  ! 

entire  signal-processing  system,  and  in  Fig.  4.5  a set  of  sample  waveforms 
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Figure  4.3.  Sketch  of  sample  waveforms  at  the  output  terminals 
of  the  high-pass  filter,  rectifier,  and  low-pass 
filter. 
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Figure  4.5.  Sample  waveforms  obtained 
at  various  points  (A,B,C,  etc.)  in 
signal  processing  chain  of  Figure  4.4. 
Upper  trace  is  current  signal,  lower 
trace  is  ^Vi^ner  signal,  excitation 
frequency  = 15  khz,  artificially- 
induced  modulation  at  = 500  hz  to 
simulate  fluctuations  in  conductivity. 
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obtained  at  each  station  in  the  signal-processing  chain,  using  a dummy 
probe  with  artificially  induced  fluctuations  in  the  i and  6V  signals. 

4.3  Procedure  for  Calibrating  the  Electronics,  Using  a Dummy  Probe 

A simplified  diagram  of  the  signal-processing  system  of  Fig.  4.4  is 
shown  in  Fig.  4.6,  in  which  the  differential  amplifier,  high-pass  filter, 
rectifier,  low-pass  filter,  and  variable-gain  amplifier  for  each  of  the 
numerator  and  denominator  chains  have  been  "lumped"  into  a single  compo- 
nent having  some  net  gain  factor,  which  equals  the  product  of  the 

individual  gain  factors  of  its  constituent  components.  In  addition.  Fig. 
4.6  shows  switches  which  allow  the  operator  to  switch  the  current  supplied 
by  the  probe  power  supply  back  and  forth  from  the  "real"  four-pin  probe 
(which  may  be  immersed  in  either  the  plasma  or  in  the  calibrating  electro- 
lyte solution)  to  a set  of  three  fixed  resistors  (R2,R^,  and  ^^^2^  which 
constitute  what  we  refer  to  as  a "dummy  probe". 

The  operator  adjusts  all  the  gain  factors  of  the  electronic  components 
such  that,  when  the  switches  connect,  alternately,  the  real  and  the  dummy 
probes  into  the  circuit,  all  the  signals  remain  within  their  allowable 
ranges.  Once  this  is  done,  the  operator  switches  to  the  dummy  probe  and 
reads  the  corresponding  constant  divider  output,  accurate 

voltmeter.  Then  he  switches  over  to  the  real  probe  and  reads  (or  records 
on  tape  or  a strlpchart)  the  corresponding  (possibly  time-varying)  divider 
output,  Qj-eal'  these  two  measurements,  plus  a knowledge  of  the  probe 

constant  k^  and  the  value  of  the  resistor  ^2.2’  conductivity  of  the 

fluid  in  which  the  real  probe  is  immersed  can  be  determined,  as  explained 


below. 
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from  which  it  follows  that 


40 


Numerator 


processing  electronics 


^net.denom 

B 

“net,num 


*^dunnny^l2 


(4.3) 


From  Eqn.  (1.5),  a = ^^isV)  ‘ after  suitable  signal  processing, 

as  discussed  above,  the  i-6V  characteristics  are  linear  and  pass  through 

the  origin  of  coordinates,  so  we  have  that  O = k (I  /6V  ),  which 

° p max  max  ’ 

may  be  written  in  terms  of  the  notation  of  Fig.  4.6  as  follows: 


I V k V /b 

. max  _ . shunt, max  _ p N,real  °net,num 

p 6V  " p R ,.6V  ° R u VT  ,/g  " j 

max  sh  max  sh  D,real  °net,denom 


1 ®net , denom  \ 

/in  N,real\ 

\ g 1 

' ®net,nuin  ' 

\ ^ V 1 

' D,real' 

(4.4) 


Substituting  (4.3)  into  (4.4)  and  identifying  the  last  term  in  brackets  in 
(4.4)  as  the  divider  output  obtain 

k Q , (t) 


a(t)  = 


k Q ,(t) 
p real 

^12  ^dummy 


(4.5) 


which  is  the  working  equation  for  determining  the  electrical  conductivity 
with  the  four-pin  probe  and  electronics  setup  already  described. 


4.4  Procedure  for  Determining  the  Probe  Constant,  k 


r 

The  probe  constant  k in  Eqn.  (4.5)  was  determined  experimentally 
using  the  Identical  electronics  setup  shown  in  Figs.  4.4  and  4.6.  The 
only  differences  in  procedure  were  that  (a)  the  probe  being  calibrated 
was  immersed  in  a carefully  prepared  aqueous  KOI  solution  whose  conduc- 
tivity, f®  accurately  known  from  handbooks,  and  (b)  the  excitation 

frequency  and  cutoff  frequencies  of  the  filters  were  different  when  using 
the  probe  in  the  electrolyte  solution  than  when  using  the  probe  in  the 
plasma.  For  measurements  in  the  electrolyte,  excitation  frequencies  around 

1 kHz  were  found  to  be  most  suitable,  and  measurements  of  k were  in  fact 

P 

taken  over  a range  of  frequencies  (typically  100  Hz  to  10  kHz).  The  high- 
pass  filter  cutoff  was  lowered  to  around  10  Hz  for  these  measurements. 

For  each  excitation  frequency,  the  following  quantities  were  measured: 
the  KCS,  solution  temperature  (which  was  needed  to  accurately  determine 
the  conductivity  from  tabulated  tata) , the  divider  output 


f 


when  the  dummy  probe  was  switched  into  the  circuit,  and  the  divider  out- 
put Qj-eal  KCi  obtained  when  the  real  probe  being  calibrated  is  immersed 
in  the  standard  KC£  solution.  Substituting  these  quantities  into  the 
working  equation  (4.5)  and  solving  for  yields 


= R 


12  KC£  Q 


Mummy 


(4.6) 


real.KCS. 


It  was  found  that. 


kp  as  determined  by  (4.6)  had  a fairly 


which  is  the  working  relation  used  to  determine  k_^ 
unless  great  care  was  taken, 
large  frequency  dependence,  related  to  polarization  effects  in  the  elec- 
trolyte. Appendix  A presents  in  more  detail  the  procedures  used  to  over- 
come these  problems  and  obtain  reasonably  accurate  values  of  k^.  The 
predominantly  "ratiometric"  character  of  the  whole  signal-processing  and 
experimental  procedure  is  illuminated  by  substituting  Eqn.  (4.6)  into 
Eqn.  (4.5)  to  obtain  the  result 


0 , (t)  = 

plasma 


^ ^^real , plasma  '^^real  ,KC)l^ 
KCX,  plasma ,KCil^ 


(4.7) 


Chapter  5 

INVESTIGATIONS  OF  FLOW  EFFECTS  ON  INDICATED  CONDUCTIVITY' 


5.1  Description  of  Experiments 

Investigations  of  flow  effects  on  the  time-averaged  conductivity  in- 
dicated by  the  four-pin  probes  were  carried  out  using  two  probes  simulta- 
neously, a fixed  probe  and  a mobile  probe,  both  operated  with  identical 
but  separate  electronics.  As  explained  in  Section  4.2,  the  outputs  of 
the  analog  dividers  for  each  probe  were  sent  through  final  low-pass  fil- 
ters (whose  cutoff  frequencies  were  usually  set  at  around  2 Hz)  to  accom- 
plish the  time  averaging,  and  then  to  separate  channels  of  an  eight- 
channel  analog  strlpchart  recorder.  The  fixed  probe,  located  in  the  low- 
velocity  region  of  the  plenum  chamber,  was  employed  mainly  because  many 
prior  experiments  had  shown  that  slow,  often  erratic,  and  often  fairly 
large  changes  in  the  conductivity  of  the  seeded,  arc-heated  argon  flow  did 
occur  even  for  nominally  constant  flow  conditions.  The  purpose  of  the 
fixed  probe  was  thus  to  monitor  these  changes  in  the  "'background"  conduc- 
tivity which  were  not  due  to  flow  effects,  so  as  to  avoid  erroneously 
interpreting  variations  of  indicated  conductivity  on  the  mobile  probe  as 
being  due  to  flow  effects.  Thus,  the  fixed-probe  output  served  as  a 
"normalization  factor",  and  in  all  the  final  results  of  these  experiments 
we  plot  the  quotient,  %obile^^f ixed ’ indicated  conductivities 

from  the  two  probes  at  the  same  instant  in  time. 

The  procedure  used  to  study  the  effects  of  flow  velocity  on  the  indi- 
cated conductivity  translate  the  mobile  probe  axially  in 

and  out  along  the  centerline  of  the  converging  nozzle,  thereby  varying  the 
flow  velocity  at  the  probe  tip.  The  axial  position,  z,  of  the  probe 
(from  which  the  velocity  was  later  deduced  as  described  in  Appendix  C)  was 
monitored  continuously  by  a position  transducer  whose  output  was  recorded 
on  the  eight-channel  strlpchart. 

Similarly,  to  study  the  effects  of  the  orientation  of  the  pins  and  to 
measure  radial  profiles  of  conductivity  in  the  flow,  the  probe  mount  was 
slowly  rotated  through  the  angle  0,  or  translated  radially  through  radial 
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distance  r normal  to  the  flow.  The  variables  r and  6 were  monitored 


with  position  transducers  and  similarly  recorded  on  the  same  eight-channel 
strlpchart. 

In  addition  to  the  indicated  conductivities  of  the  two  probes  and  the 
position  variables  z,  r,  and  9,  several  other  signals  were  simulta- 
neously recorded  on  the  strlpchart.  These  signals  included  the  gas  tem- 
perature (actually  the  junction  emf  of  a tungsten  thermocouple  in  the 
plenum),  the  arcjet  current  and  voltage  drop,  and  a voltage  proportional 
to  the  rate  of  Injection  of  seed  (NaK)  into  the  plenum. 


5.2  Effects  of  Flow  Velocity  on  Indicated  Conductivity 

Figures  5.1  and  5.2  are  reproductions  of  portions  of  actual  strip- 
chart  records  obtained  in  the  experiments  just  described.  These  figures 
show  the  output  of  the  axial  position  transducer  and  of  the  fixed  and  mo- 
bile probes,  and  illustrate  that,  while  the  fixed  probe  output  remains 
fairly  constant,  the  mobile  probe  output  decreases  as  the  velocity  increases. 
In  other  words,  these  strlpchart  records  already  indicate  the  existence  of 
a substantial  velocity  effect  on  the  indicated  conductivity.  Figs.  5.1  and 
5.2  show  two  different  "modes"  of  data-taking:  in  Fig.  5.1  the  mobile  probe 
was  continuously  withdrawn  and  inserted  (hence  the  name  "continuous"  mode) , 
whereas  in  Fig.  5.2  the  mobile  probe  was  set  at  several  discrete  axial  sta- 
tions and  data  were  recorded  for  at  least  30  seconds  at  each  station  ("dis- 
crete" mode) . Because  the  probes  may  have  some  short  "transient  response 
time"  after  a velocity  change,  during  which  the  probe  output  is  in  error, 
the  data  obtained  in  the  "discrete"  mode  may  possibly  be  more  accurate 
than  the  "continuous"  mode  data. 

After  reducing  these  and  other  similar  strlpchart  data,  employing 

formula  (4.5)  to  obtain  ’^fixed’  plotting  the  quotient 

of  these  results  vs.  u/u  , we  obtain  the  curve  shown  in  Fig.  5.3,  in 

max 

which  only  "discrete"  mode  data  points  were  used  in  sketching  in  the  con- 
tinuous curve.  As  can  be  seen,  the  indicated  conductivity,  '^niobile’ 
reduced  by  around  50%  as  the  velocity  Increases  over  the  range  of  veloci- 
ties achievable  in  these  experiments  (namely  from  about  10  m/sec  to  about 
100  m/sec).  This  result,  while  obtained  using  a particular  size  four-pin 
probe  in  the  "end-on"  flow  configuration,  nevertheless  indicates  that  there 
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respectively,  to  the  indicated  electrical  conductivities 
‘^fixed  %obile)*  Lower  trace  shows  the  axial 
posirion  of  the  mobile  probe  being  set  at  several 
discrete  positions  in  the  nozzle. 


Figure  5.3.  Ratio  of  the  indicated  conductivity  of  a mobile  4-pin  probe  to  that  of  a fixed 
position  probe,  as  a function  of  flow  velocity  u at  the  mobile  probe  tip. 

The  continuous  curve  was  sketched  using  only  the  "discrete"  mode  data.  See 
Figure  5.4  for  mobile  probe  dimensions. 


is  a significant  "flow  effect"  (vs.  a velocity  dependence  in  the  indicated 
conductivity)  which  must  be  taken  into  account  when  using  four-pin  probes. 

In  interpreting  the  results  shown  in  Fig.  5.3  as  evidence  of  a "flow 

effect",  we  have  assumed  that  the  local  conductivity  remains  constant  in 

the  core  flow  of  the  plasma  between  the  upstream  and  downstream  ends  of 

the  nozzle.  In  Appendix  F we  discuss  in  more  detail  this  assumption. 

While  there  are  mechanisms  (such  as  gasdynamics  cooling  and  radiation 

escape)  which  would  tend  to  lower  a,  the  conclusion  of  Appendix  F is 

that  it  is  likely  that  the  electron  density  (and  hence  o)  was  partially 

frozen  at  its  upstream  value  as  the  gas  flowed  down  the  nozzle  in  these 

experiments.  Also  included  in  Appendix  F are  estimates  of  the  conductivity 

drops  which  would  occur  from  radiation  escape  and  gasdynamic  cooling,  even 

in  the  absence  of  frozen  flow,  and  these  results  all  predict  less  than  a 

14%  conductivity  drop,  well  below  the  indicated  50%  drop  that  was  measured. 

Still  another  mechanism  leading  to  a real  conductivity  drop  would  be  a 

cooling  of  the  gas  through  heat  transfer  to  the  nozzle  walls.  However, 

direct  measurements  of  the  radial  profiles  of  conductivity,  presented  below 

in  Section  5.4  (cf.  Fig.  5.8)  indicate  that  there  is  a large  core  region  in 

the  flow  for  which  3o/3r  = 0.  This  is  fairly  strong  evidence  in  support  of 

the  claim  that  9T/3r  = 0 in  the  core  flow  region,  and  hence  that  this  region 

is  not  cooled  by  heat  conduction.  Thus  the  mechanisms  of  gasdynamic  cooling, 

radiation  escape,  and  heat  conduction  are  not  sufficient  to  explain  the 

observed  decrease  in  o.  . 

indic 

A convenient  way  of  interpreting  the  results  shown  in  Fig.  5.3,  and  one 

which  allows  for  a quantitative  correction  factor  to  be  introduced  for  the 

effect  of  flow  velocity  on  four-pin  probe  readings,  is  to  regard  the  probe 

constant  k as  a velocity-dependent  quantity.  According  to  the  working 
P 

equation  (4.5),  we  have  that 


k = (R,„Q,  ) ~ 

p 12  dummy  Q 


(5.1) 


and,  since  R, „ , Q , , and  o are  constants  in  these  experiments,  the 

12  dummy’ 

decrease  in  divider  output,  Q,  as  velocity  increases  can  be  thought  of 
as  an  increase  in  the  probe  constant.  (a  may  in  reality  vary  by  a few 
percent  along  the  nozzle  axis,  but  for  simplicity  we  assume  it  to  be  essen- 
tia-ly  constant.)  The  quantity  of  interest,  then,  for  a given  probe,  is 


i 
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the  ratio  of  the  values  of  the  probe  constant  at  zero  velocity  (a  quantity 
measurable  to  acceptable  accuracy  In  a salt  water  solution)  to  that  at 


some  non-zero  velocity,  u.  For  convenience  we  Introduce  the  notation 


f(u)  = 


kp(u) 

k (u=0) 
P 


(5.2) 


for  this  ratio,  and  we  shall  show  how  f(u)  can  be  determined  from  the 

experimental  curve  of  a,  (u)/a.  ,(u  , )vs.  u/u 

Indlc, mobile  Indic, fixed  plenum  max 

shown  In  Fig.  5.3. 

By  writing  Eqn.  (5.1)  for  the  mobile  probe  at  some  arbitrary  velocity 
u,  and  dividing  by  the  same  equation  written  for  the  fixed  probe  at  the 
particular  non-zero  but  constant  velocity  ^pj^enum  present  In 

the  plenum  at  the  location  of  the  fixed  probe,  we  obtain  the  result 


k , T (u) 
, mob lie 


k j(u  1 

p, fixed  plenum 


’ (, 


^oblle^*^^ 


^flxed^%lenum 


'12  ,moblle^dummy,moblley  '^2  , f Ixed^dummy , f Ixed 
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Multiplying  both  sides  of  this  equation  by  the  constant  ratio 

k,.  j(u=0)/k  . (u=0)  leads  to  the  result 

fixed  mobile 


/ ^p,moblle^^^  \ 

\k  ...  (u=0)  / 
' p, mob  lie 


/ *^p , f Ixed  ^ ^p lenum^  \ 
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^moblle^^^ 


f Ixed^^plenum^ 
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^12, 
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Mobile 
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dummy, mob lie 


:]/P 


,(u=0)  0,.  ,(u  , )■ 

),  fixed fixed  plenum 


12, fixed 


'dummy , f Ixed 


In  which  the  terms  In  square  brackets  are  just  '^j^ndlc’  calculated 
from  experimental  data  using  Eqn.  (4.5).  Hence  we  arrive  at  the  result 


f , (u) 

mobile 


^flxed^^plenum^  ^lndlc,moblle^'^^  ^ '^lndlc,flxed^'^plenum^ 


(5.3) 


In  which  the  right-hand  side  Is  the  numerical  Inverse  of  the  curve  plotted 
In  Fig.  5.3. 
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We  still  do  not  know  the  value  of  the  constant  f,.  j(u  , ), 

fixed  plenum 

but  it  can  be  found  if  we  permit  ourselves  to  extrapolate  the  curve  in 

Fig.  5.3  to  u=0.  At  u=0,  we  have  that  f , (u=0)  = 

® ’ mobile 

S,moblle^“=°^/S,moblle^^=°^  = 


^fixed^^plenum^ 


^moblle^^=Q^ 

^f ixed^^plenum^ 


/vertical  intercept  of \ 
\Fig.  5.3  at  u = 0 ' 


which  has  in  this  case  the  value  f,.  ,(u  . ) = 1.4  ± 0.1,  with  the 

fixed  plenum  ’ 

uncertainty  being  estimated  from  the  uncertainty  possible  in  extrapolating 
Fig.  5.3  to  u = 0.  With  this  result  we  have  from  (5.3)  that 


^mobile^“>  " - 


„)) 


indie, mobile  / indie, fixed  plenum 
= (1.4  ± 0.1) (Inverse  of  the  curve  of  Fig.  5.3)  . (5.4) 

Using  (5.4)  together  with  the  experimental  curve  in  Fig.  5.3,  the 
corresponding  curve  of  the  probe  constant  ratio  ^mobile^'^^  calculated 

as  a function  of  velocity,  and  is  plotted  in  Fig.  5.4.  This  curve  allows 
a correction  for  velocity  effects  to  be  made.  The  conductivity  indicated 
by  the  four-pin  probe  is  to  be  calculated  from  the  formula 


k (u=0) 


R 


f(u) 


12 


■‘dummy 


(5.5) 


or  equivalently 


a = 


kp(u=0)  f(u)  ^ 


(5.6) 


max 


in  which  kp(u=0)  is  the  value  found  using  stationary  electrolytes  and 
f(u)  is  read  from  Fig.  5.4. 

It  should  be  noted  that  the  curve  for  f(u)  shown  in  Fig.  5.4  was 
obtained  using  only  one  particular  probe,  whose  dimensions  are  shown  in 
Fig.  5.4,  and  this  probe  was  employed  in  the  "end-on"  orientation  with  re- 
spect to  the  flow.  Thus,  while  the  present  work  shows  that  it  is  possible 
to  take  velocity  effects  into  account  for  any  given  probe  by  performing  ex- 
periments similar  to  the  ones  described  in  this  report,  it  is  not  known  to 
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what  extent  the  correction  factor  f(u)  may  vary  with  other  probe  geom- 
etries or  flow  conditions.  Furthermore,  since  the  correction  factor 
f(u)  is  relatively  large  (varying  between  unity  and  nearly  3 in  the  pres- 
ent work) , it  would  appear  quite  important  to  perform  such  detailed 
"velocity-effect  calibration"  measurements  on  any  four-pin  probe  used  in 
MHD  conditions,  if  accurate  absolute  values  of  a are  desired.  Failure 
to  take  velocity  effects  into  account  at  all  could  result,  according  to 
the  data  of  this  investigation,  in  indicated  conductivities  which  are  as 
much  as  a factor  of  two  or  three  lower  than  the  actual  conductivity. 

5.3  Effects  of  Pin  Orientation  on  Indicated  Conductivity 

A separate  experiment  was  carried  out  to  investigate  the  possible 
effects  of  the  orientation  of  the  probe  pins  relative  to  the  flow  direc- 
tion, and  several  difficulties  were  encountered  which  resulted  in  consid- 
erably less  data,  and  somewhat  poorer  quality  data,  than  in  the  velocity 
effects  experiments.  Also,  in  order  to  disturb  the  gas  flow  field  as 
little  as  possible,  the  particular  mobile  probe  used  for  these  experiments 
was  made  much  smaller  than  any  other  previously  used  probes;  the  double- 
bore alumina  tubing  used  had  a 1/16-inch  outer  diameter,  the  tungsten  wire 
diameter  was  0.011  inches,  and  dimension  "a"  (see  Fig.  1.1)  was  0.022  inch 
(so  that  the  separation  of  the  metal  surfaces  of  the  adjacent  outer  and 
inner  pins  was  only  0.011  inches,  or  0.28  mm).  The  thin  ceramic  tubing 
was  somewhat  susceptible  to  warping  from  thermal  gradients,  and  during  the 
test  the  two  adjacent  ceramic  tubes  separated  by  roughly  one  diameter, 
thereby  changing  the  effective  probe  constant  (i.e.,  distance  b in- 
creased — see  Eqn.  (1.6)  and  Fig.  1.1).  This  meant  that  the  quantity 
^mobile^*^f ixed  differ  considerably  from  unity  due  to  the  uncertainty 

in  kp  Eqn.  (4.5)).  To  make  the  results  more  meaningful, 

therefore,  the  measurements  were  normalized  by  the  value  of  o , /o..  , 

’ ■'  mobile  fixed 

at  0 = 0°,  and  the  relative  results  05tO°^^°m^°f  ^6=0°  ^ 

plotted. 

A typical  stripchart  record  is  reproduced  in  Fig.  5.5,  showing  three 
relatively  fast  (~  30  second  duration)  sweeps  of  6 through  ±45“  and  the 
corresponding  drops  in  ® increases.  After  reducing  these 

data,  plus  data  from  another  much  slower  scan  of  0 to  30“  and  back  to  0“, 
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we  obtained  the  plots  shown  in  Fig.  5.6.  As  can  be  seen  in  these  plots 
(in  which  the  time  sequencing  of  the  points  plotted  is  indicated  by  the 
dashed  lines  and  arrows) , the  indicated  conductivity  curve  does  not  re- 
trace exactly  its  original  path  as  6 is  swept  out  to  some  finite  value 
and  then  back  to  zero.  There  also  appears  to  be  a slight  asymmetry  about 
the  line  0=0®.  Despite  these  unexplained  features,  the  data  seem  to  be 
sufficiently  repeatable  and  of  sufficient  quality  to  support  the  conclusion 
that  as  the  orientation  of  the  pins  is  changed  from  0=0°  ("end-on" 
orientation  with  the  pins  parallel  to  the  flow)  to  0 = 45°,  the  indi- 
cated conductivity  drops  by  20%  to  25%,  and  this  decrease  appears  to  be 
linear  in  0,  at  least  to  a first  approximation. 

5.4  Spatial  Resolution  and  Effects  of  Nearby  Conductors 

We  mentioned  in  Section  1.3  that  the  spatial  resolution  of  a four-pin 
probe  depends  on  the  size  of  the  probe.  For  purposes  of  defining  more  pre- 
cisely what  we  mean  by  "spatial  resolution",  let  us  suppose  that  the  probe 
tip  is  surrounded  by  a sphere  of  uniform-conductivity  plasma  of  radius  r, 
and  that  beyond  this  radius  the  conductivity  may  be  different,  or  there 
may  be  conductors  or  other  disturbances.  Then  if  r is  the  largest  radius 
for  which  such  disturbances  produce  measurable  changes  in  the  conductivity 
indicated  by  the  probe,  we  say  that  the  probe  has  approximately  the  spatial 
resolution  r.  For  a four-pin  probe  it  seems  plausible  that  the  spatial 
resolution  would  be  roughly  the  dimensions  of  the  region  through  which  most 
of  the  current  between  the  outer  pins  flows.  For  probes  constructed  such 
that  the  dimension  i (see  Fig.  1.1)  is  roughly  equal  to  the  outer  pin 
separation  (a  + b in  Fig.  1.1),  we  might  anticipate  that  the  spatial  reso- 
lution would  be  approximately  equal  to  this  outer-pin  separation,  or  per- 
haps up  to  twice  this  value. 

Experimentally,  the  spatial  resolution  was  investigated  by  traversing 
the  probe  radially  across  the  flow  in  an  attempt  to  measure  radial  conduc- 
tivity profiles  in  the  nozzle.  Fig.  5.7  shows  portions  of  the  stripchart 
record  obtained  as  the  small  probe  described  in  Section  5.3  was  traversed 
several  times  across  the  flow.  There  was  some  problem  with  thermal  warping 
of  the  probe,  so  that  there  is  considerable  uncertainty  in  the  actual  posi- 
tion of  the  probe  tips  in  these  traverses,  particularly  as  the  probe 
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Figure  5. 


-*0® -30* -20*  -10®  0 10®  20®  30*  40® 

$ 

Results  of  9-scans  of  6/1/77  showing  variation  of 
indicated  conductivity  with  angle  of  orientation,  0. 
The  quantity  plotted  on  the  ordinate  is  the  ratio 


^mobile  ^mobile 

^fixed  °fixed 


(e=o). 


The  top  curve  shows  data  from  a slow,  continuous 
scan  of  6,  while  the  scans  in  the  lower  three  curves 
were  carried  out  more  rapidly. 
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approached  the  wall  where  temperature  gradients  were  higher.  Nevertheless, 
the  stripchart  shows  a relatively  flat  profile,  with  the  indicated  conduc- 
tivity increasing  very  near  the  wall. 

This  effect  is  seen  more  clearly  in  Fig.  5.8,  in  which  the  quantity 
a(rj^O) /a(r=o)  has  been  calculated  from  one  of  the  scans  in  Fig.  5.7  and 
plotted  vs.  r.  At  distances  of  about  0.15  1 0.05  inch  from  the  nozzle 
wall,  the  very  flat  profile  of  turns  sharply  upward,  in  spite  of 

the  fact  that  the  actual  local  conductivity  near  the  cooler  nozzle  wall 
must  be  decreasing.  This  effect  is  thought  to  be  due  to  a conducting  sur- 
face layer  on  the  nozzle,  and  is  consistent  with  the  observations  of  Brown 
(Ref.  [6],  p.  176),  who  found  that  nearby  conducting  walls  affected  the 
probe  current  pattern  in  such  a way  as  to  make  the  indicated  conductivity 
increase.  In  the  present  experiments,  the  conducting  layer  was  probably 
foxrmed  by  tungsten  carried  downstream  from  the  arc jet  cathode  and  reacting 
with  the  boron  nitride  nozzle  (such  layers  have  been  found  to  be  conduct- 
ing in  previous  work  at  Stanford).  Another  possible  source  for  a conduct- 
ing layer  would  be  liquid  metal  NaK  (although  the  nozzle  surface  tempera- 
ture was  supposedly  kept  high  enough  to  avoid  NaK  condensation,  as 
mentioned  in  Section  2.2).  Some  conducting  compounds  formed  by  seed  de- 
posits from  the  several  preceding  tests,  which  all  employed  the  same  nozzle, 
are  another  possibility. 

If  we  take  as  the  spatial  resolution,  r,  the  distance  between  the 
probe  and  the  wall  at  which  the  effects  of  the  conducting  wall  begin  to 
alter  (viz.,  the  figure  0.15  ± .05  inch  cited  previously),  and 

divide  by  the  separation  between  the  outer  probe  pins  (which  was  0.125 
± .03  inch),  we  obtain  the  result 

spatial  resolution,  r ^ 12+05 

outer  pin  separation  . - . , 

which  is  roughly  in  agreement  with  our  earlier  expectation  that  the  spatial 
resolution  of  the  probe  is  approximately  equal  to  the  outer  pin  separation. 

Some  additional  modeling  of  the  possible  effects  of  nearby  conductors 
and  insulators  on  the  Indicated  reading  of  four-pin  probes  was  done  using 
two-dimensional  resistance  paper  ("Teledeltos  paper").  The  base  case  was  a 
34  inch  x 34  inch  square  with  insulating  boundaries  and  electrodes  painted 
on  with  conducting  paint  (see  Fig.  5.9).  Successively  smaller  circles  were 


Resistance  Paper 


Insulating  or 

Conducting 

Boundary 


(b)  Test  cases.  In  which  r was  varied. 


Figure  5.9.  Geometry  of  modelling  of  the  effects  of  nearby 
conductors  and  insulators  with  two-dimensional 
resistance  paper. 


then  cut  out  of  this  square,  first  with  a conducting  boundary  painted  all 
around,  and  then  with  this  conducting  ring  cut  away  to  leave  an  Insulated 
boundary.  In  each  case  the  indicated  conductivity  (actually  the  slope  of 
the  i-6V  plot)  was  measured.  The  results  of  this  procedure  are  shown 
in  Fig.  5.10.  For  this  particular  model  situation  we  can  conclude  that 
the  error  in  the  indicated  conductivity  is  less  than  5%  provided  (a)  r/2W 
>2  in  the  case  of  an  insulating  boundary,  and  (b)  r/2W  >4  in  the  case 
of  a conducting  boundary.  These  results  are  of  course  dependent  on  the 
particular  geometrical  configuration  employed  for  modeling  (i.e.,  two- 
dimensional,  four-pin  probes  centered  in  circular  domains),  but  neverthe- 
less they  support  the  conclusion  that  the  spatial  resolution  of  a four-pin 
probe  is  of  the  order  of  a few  times  the  outer  pin  separation. 
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Figure  5.10.  Results  of  two-dimensional  resistance  paper  study 
of  the  effects  of  nearby  conducting  and  insulating 
boundaries  on  indicated  conductivity  from  4-pin 
probes. 
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Chapter  6 

SAMPLE  TIME-RESOLVED  CONDUCTIVITY  MEASUREMENTS 

6.1  Introduction 

The  purpose  of  this  chapter  Is  to  present  some  sample  data  taken  when 
the  four-pin  probe  was  operated  In  the  time-resolved  mode,  as  discussed  In 
Chapter  4.  The  data  are  not  Interpreted  In  the  present  work,  but  are  pre- 
sented to  Illustrate  the  capabilities  and  type  of  time-resolved  measure- 
ments possible  with  this  probe  system.  Some  Indication  that  there  were 
large  conductivity  fluctuations  present  In  the  flow  has  already  been  pre- 
sented In  Fig.  4.2,  which  shows  marked  changes  In  the  slope  of  the  1-6V 
characteristics  (and  hence  In  a)  during  the  l/60th  second  photographic 
exposure  time,  but  the  data  to  follow  In  the  next  section  show  the  time 
development  of  such  fluctuations  directly,  using  time-base  swept  oscillo- 
scope photographs. 

6.2  Sample  Oscilloscope  Photographs 

Figure  6.1a  shows.  In  the  upper  two  traces,  the  waveforms  of  the  am- 
plified 1 and  6V  signals  as  they  emerge  from  the  high-pass  filters 
(l.e.,  at  point  C In  Fig.  4.4).  The  lower  trace  shows  the  voltage  be- 
tween the  two  outer  pins,  '^'^outer'  Several  features  of  this  photograph 
are  worth  commenting  upon:  (1)  Although  the  applied  voltage  at  the  outer 
pins  Is  sinusoidal  (or  nearly  so) , the  1 and  6V  signals  are  clearly 
distorted  Into  a nonslnusoldal,  almost  triangular  waveform.  This  effect 
Is  not  attributable  to  the  signal-processing  electronics,  since  In  the  cor- 
responding photograph  (not  Included  here)  of  the  1 and  6V  signals  from 
the  dummy  probe  (cf.  Section  4.3),  all  the  waveforms  are  sinusoidal.  Thus 
It  Is  suspected  that  the  non-slnusoldal  1 and  6V  signals  arise  from 
some  "real"  effect  In  the  plasma.  (11)  Despite  the  difference  In  waveforms 
between  the  ^^gy^er  ^^snal  on  the  one  hand  and  the  1 and  6V  signals 
on  the  other,  the  1 signal  does  appear  to  have  the  same  (distorted)  wave- 
form as  the  6V  signal,  so  that  the  ratio  of  the  amplitudes  of  these  two 
signals  will  still  be  proportional  to  the  slope  of  the  1-6V  characteristic. 
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Figure  6. 1 . 


(C) 


Sample  Oscilloscope  Photographs,  (a)  and  (b).  Upper 
Trace:  amplified  current  signal  (2  V/div);  Middle 
Trace:  amplified  inner  pin  voltage  difference  (2  V/div); 
Lower  Trace:  outer  pin  voltage  difference  (1  V/div). 

Photo  (a):  from  mobile  probe.  Run  90,  3/15/77;  sweep  speed 
10  vjsec/div;  excitation  frequency  = 40  kHz; 

Photo  (b):  from  fixed  probe.  Run  130,  2/11/77;  sweep  speed 
100  usec/div;  excitation  frequency  = 10  kHz 
(band-oass  filter  used  in  place  of  high-pass 
filter  shown  in  Figure  4.4); 

Photo  (c):  from  fixed  probe.  Run  410,  3/15/77.  Upper  Trace 
(with  respect  to  upper  most  horizontal  basel ine, 
5-V/div):  probe  current  signal  at  input  to 
analog  divider.  Middle  Trace  (with  respect  to 
middle  baseline,  5 V/div) : fnner  pin  voltage 
difference  signal  at  input  to  divider.  Lower 
Trace  (with  respect  to  lowest  baseline,  5 V/div): 
analog  divider  output,  proportional  to  a(t). 

Sweep  speed  5 msec/d iv  in  photo  (c). 
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a multiplicative  constant,  is  equal  to  the  plasma's  conductivity  as  a 
function  of  time.  This  photo  clearly  shows  the  presence  of  very  large 
fluctuations  in  O in  the  plenum  chamber,  where  the  probe  was  located. 

It  also  gives  some  idea  of  the  typical  time  scales  of  the  fluctuations, 
and  demonstrates  that  the  analog  divider  is  working  properly. 

6.3  Sample  Statistical  Data  on  Conductivity  Fluctuations 

In  this  section  we  Include  some  statistical  properties  of  the  elec- 
trical conductivity  fluctuations  which  were  present  in  the  seeded  argon 
flow  employed  in  this  research.  The  procedure  used  to  obtain  this  quanti- 
tative statistical  information  was  first  to  make  analog  tape  recordings 
of  the  processed  output  signal  (proportional  to  o(t))  from  each  of  the 
four-pin  probes  (as  described  in  detail  in  Section  4.2),  and  then,  subse- 
quent to  the  experiment,  to  perform  numerical  computations  on  these  recor- 
ded data.  The  computations  were  accomplished  by  first  sending  the  recorded 
probe  output  (proportional  to  a(t))  into  a Preston  GMAD-1  Analog-to- 
Digital  Conversion  System,  and  then  using  this  digital  data  as  input  for 
various  statistical  programs  executed  by  an  HP-2100  minicomputer.  The 
computed  results  included  the  following  quantities:  The  mean  conductivity, 
a;  the  rms  value  of  the  conductivity  fluctuations,  “ |^E(a(t)-o)^J  S 

the  spectral  density  function  (or  "frequency  spectrum");  the  probability 
density  function  (all  of  these  for  recordings  from  individual  probes  at 
specific  points  in  the  flow);  and  the  cross-correlation  coefficient  between 
the  outputs  of  two  probes  in  different  locations.  In  each  case  the  propor- 
tionality constant  between  the  conductivity  Cj(t)  and  the  recorded  probe 
output  voltage  Q(t)  was  corrected  for  the  effect  of  velocity  according 
to  the  procedure  described  in  the  preceding  chapter  (cf.  Eqn.  (5.5)). 

This  procedure  assumes  that  the  velocity-effect  correction  factor,  deter- 
mined from  time-averaged  data,  remains  valid  for  fast  time-resolved  mea- 
surements also. 

Figure  6.2  is  a plot  of  the  ratio  of  io  o vs.  velocity. 

This  figure  Indicates  that  the  flow  velocity  (or  alternatively  the  axial 
position  in  the  converging  nozzle)  does  not  appreciably  affect  the  rela- 
tive conductivity  fluctuation  level.  In  fact,  the  absolute  value  of  o 

rms 

itself  is  quite  constant  for  these  data,  since  the  corrected  mean 
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as  we  require  for  the  signal- processing  system  to  yield  an  output  which 
is  truly  proportional  to  a.  (ill)  There  is  some  phase  shift  between 
6Vouter  ^ signals.  This  is  due  to  the  signal- 

processing  components  which  shift  the  phase  of  i and  6V  by  the  same 
amount  relative  to  that  of  6V  ; the  same  phase  shift  occurs  when  the 
dummy  probe  is  used.  This  phase  shift  causes  no  difficulties  so  long  as 
it  is  the  same  for  both  the  i and  6V  signals,  as  indeed  seems  to  be 
the  case  in  Fig.  6.1a.  \ 

Figure  6.1b  displays  similar  signals  as  Fig.  6.1a,  but  from  a differ- 
ent test  and  using  a different  probe  and  slightly  modified  signal- 
processing components.  In  particular,  in  place  of  the  high-pass  filter 
shown  in  Fig.  4.4,  a "high-Q"  band-pass  filter,  tuned  to  the  source  exci- 
tation frequency,  was  used.  This  configuration  worked  well,  but  gave  the 
somewhat  misleading  impression  that  the  i and  6V  signals  were  sinus- 
oidal, as  was  the  excitation  waveform.  In  fact,  as  seen  in  Fig.  6.1a, 
this  is  not  the  case.  The  reason  for  the  "sinusoidal-looking"  waveforms 
in  Fig.  6.1b  is  that  the  band-pass  filters  attenuate  any  high-frequency 
Fourier  components  of  the  input  waveform,  and  it  is  just  these  higher  fre- 
quency components  which  give  the  waveforms  of  Fig.  6.1a  their  "sawtooth" 
character. 

The  main  reason  for  including  Fig.  6.1b  in  this  report  is  to  illus- 
trate the  "slow"  changes  in  the  amplitudes  (or  modulation)  of  the  higher- 
frequency  "carrier"  signals,  i and  6V.  This  effect  is  clearly  seen  in 
the  waveforms  of  Fig.  6.1b;  the  modulation  is  caused  by  relatively  much 
slower  fluctuations  in  the  plasma  conductivity  which  occurred  during  the 
one  millisecond  in  which  the  oscilloscope  trace  swept  across  the  screen. 

The  same  effect  is  to  be  expected  in  Fig.  6.1a,  but  does  not  show  up  so 
clearly  due  to  the  much  faster  sweep  speed  employed. 

The  sample  traces  shown  in  Fig.  6.1c  are  Included  to  convey  a visual 
impression  of  typical  fluctuations  in  the  rectified  probe  current  i,  inner 
pin  potential  difference  6V,  and  conductivity.  The  upper  trace  is  the 
"completely  processed"  current  signal  (appearing  at  point  "F"  in  Fig.  4.4) 
which  is  sent  to  the  numerator  input  terminal  of  the  analog  divider;  the 
middle  trace  is  the  correspondingly  processed  6V  signal  which  goes  to  the 
denominator  input.  The  lower  trace  is  the  divider  output,  which,  aside  from 


conductivity  a = 11.5  mho/m  only  varied  by  less  than  10%  among  the 
points  plotted  In  Fig.  6.2.  One  possible  explanation  which  would  be  con- 
sistent with  these  results  Is  that  the  regions  of  nonuniformity  In  the 
plasma  are  simply  convected  downstream  with  little  diffusion  or  mixing 
during  the  approximately  3 millisecond  residence  time  In  the  nozzle.  If 
appreciable  mixing  of  the  nonuniformities  did  occur,  making  the  plasma 
more  uniform  near  the  end  of  the  nozzle  (l.e..  In  the  hlgh-veloclty  region), 

one  would  expect  the  ratio  0 /a  to  decrease.  The  fact  that  a , 

rms  rms 

computed  using  the  velocity-effect  correction  factor  f(u),  does  not  In- 
crease as  the  gas  flows  down  the  nozzle  tends  to  support  the  Interpreta- 
tion and  use  of  f(u)  and  the  assumption  that  while  f(u)  was  calculated 
from  time-averaged  data  It  Is  valid  for  time-resolved  measurements  also. 

Figure  6.3  shows  semi- log  plots  of  the  power  spectral  density  func- 
tion (see,  for  example.  Reference  [9]),  G^(f),  defined  as 

G^(f)  = 11m  \ f [o(t,f , Af)  - a]  dt  . 

Af-H)  J o 

This  quantity  represents  the  Fourier  spectrum  of  the  mean  square  conduc- 
tivity fluctuation.  It  is  clear  from  these  plots  that  the  dominant  con- 
tribution to  the  conductivity  fluctuations  occurs  at  frequencies  below 
about  200  Hz,  and  that  there  is  negligible  contribution  from  frequencies 
over  2 kHz.  The  plot  shown  for  the  case  where  the  mobile  probe  was  located 
in  the  high-velocity  region  near  the  nozzle  exit  lies  above  the  correspond- 
ing plot  for  the  low-velocity  region.  This  behavior  Is  generally  consis- 
tent with  the  interpretation  that  the  nonuniform  regions  of  fluid  are 
convected  downstream  past  the  probes,  and  that  as  the  fluid  moves  faster 
the  same  nonuniformities  give  rise  to  a higher-frequency  component  in  the 
probe  reading.  There  seems  also  to  be  an  "extra"  very  low-frequency  com- 
ponent in  the  signal  from  the  probe  when  near  the  exit.  This  additional 
low-frequency  noise  showed  up  on  the  eight-channel  stripchart  as  well, 
and  may  be  due  to  the  close  proximity  of  shear  layers  and  low-frequency 
vortex-shedding  effects  at  the  nozzle  exit,  whose  Influence  could  easily 
be  felt  in  the  subsonic  flow  at  the  location  of  the  probe,  which  was  ap- 
proximately 0.5  cm  upstream  from  the  exit  plane. 
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Figure  fi.3.  Log  of  spectral  density  function  of  conductivity  fluctuations 
vs.  freguency  from  data  of  ^121111. 


A sample  probability  density  function  for  the  electrical  conductivity 
is  shown  in  Fig.  6.4.  The  most  striking  feature  of  this  curve  is  the  ob- 
vious skewness,  which  means  that  small,  negative  deviations  of  a from 
a are  more  likely  than  positive  deviations,  but  that  large  positive  de- 
viations are  more  probable  than  large  negative  deviations.  This  behavior 
is  qualitatively  consistent  with  the  hypothesis  that  the  temperature  fluc- 
tuations in  the  plasma  have  a symmetrical  distribution  about  the  mean.  If 
Saha  equilibrium  obtains  in  the  plasma,  then  there  will  be  a highly  non- 
linear but  nevertheless  one-to-one  correspondence  between  O and  T^ 
such  that  a symmetric  temperature  probability  density  will  yield  a skewed 
conductivity  probability  density  of  the  general  shape  of  Fig.  6.4.  In 
fact,  the  temperature  probability  density  function  which  corresponds  to 
Fig.  6,4  could  easily  be  reconstructed.  If  f(a)  is  the  probability  den- 
sity for  a and  a = 8(Tg)  is  the  functional  relationship  between  O 
and  T^,  then  the  probability  density  h(T^)  for  the  temperature  is 
given  by  the  relation 


h(T^)  = f[g(Tg)]  g'(T^)  . 


The  final  statistical  calculation  which  we  shall  discuss  here  deals 

with  the  cross-correlation  coefficient  p (t)  s „ (t)/V  , where 

m f m f f m 

the  cross-correlation  function  R (t)  between  the  indicated  conductivi- 

O 

ties  a and  0^  of  the  mobile  and  fixed  probes  is  defined  as 
m f 

T 

“o,.,  = if  • 

f m •'o 


and  V and  V are  the  root  mean  square  values  of  cr,(t)  and  a (t). 
0.0  t m 

f m 

These  calculations  were  done  as  part  of  an  attempt  to  determine  experimen- 
tally the  flow  velocity  as  a function  of  axial  distance  along  the  nozzle. 

The  basic  idea  of  this  technique  was  as  follows;  With  the  fixed  and  mobile 
probes  separated  by  some  known  distance,  S,  the  output  signals  from  both 
probes  were  simultaneously  recorded  on  the  same  multi-channel  analog  tape 
for  a suitable  time  period  (e.g.,  30  seconds).  From  these  data  the  cross- 
correlation coefficient  ^ (t)  was  computed.  Assuming  that  any  fluctu- 

m f 

ations  seen  by  the  upstream  fixed  probe  are  essentially  convected  downstream 
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Figure  6.4.  Probability  density  function  for  electrical  conductivity,  mobile 
probe  in  high  velocity  region  run  102,  5/27/77,  u/u  = 0.94. 


I 

I 

1 

1 

and  seen  at  a slightly  later  time  hy  the  mobile  probe,  there  should  he 

a peak  in  the  function  some  delay  time  x^,  which  is 

equal  to  the  transit  time  of  the  nonuniformities  as  they  travel  between 

; the  two  probes.  From  several  such  recordings  and  calculations  a plot 

I could  be  made  of  the  distance  s traveled  by  a fluid  particle  as  a func- 

I tion  of  the  transit  time  required,  and  the  flow  velocity  at  any 

; location  determined  from  the  slope  ds/dx  . 

I s 

In  practice  it  did  not  prove  feasible  to  make  an  accurate  determi- 
i nation  of  velocity  in  the  manner  just  outlined,  but  the  qualitative  trends 

! were  as  expected.  A typical  cross-correlation  coefficient  plot  is  shown 

in  Fig.  6.5.  There  is  a rather  broad  peak  and  a corresponding  transit 
I time  x^  = 1.5  msec.  However,  the  broad  peak  Indicates  that  there  is  a 

j high  correlation  between  the  two  signals  over  quite  a range  of  delay  times 

I X,  or,  in  other  words,  the  physical  scale  length  of  the  nonuniformities 

must  be  large  compared  to  the  separation  between  the  probes.  Also,  tape 
recordings  were  made  at  only  a few  separations  in  the  5 cm  to  10  cm  range, 
and  the  cross-correlation  coefficient  curves  showed  less  correlation  and 
even  broader  peaks  as  the  separation  increased.  The  limited  and  somewhat 
1 uncertain  data  on  s vs.  X^  which  was  obtained  is  plotted  in  Fig.  6.6. 

While  the  delay  times  do  increase  with  separation  generally,  there  are 
three  points  showing  the  same  delay  time.  Also,  for  the  very  largest 
separation  employed,  the  peak  in  the  cross-correlation  coefficient  was  so 
broad  and  noisy  that  it  was  very  arbitrary  which  delay  time  to  use.  Hence 
it  was  not  possible  to  obtain  accurate  velocities  by  this  technique.  A 
rough  calculation  of  the  slope  of  a line  drawn  through  the  three  "best" 
points  in  Fig.  6.6  (see  the  dashed  line)  yields  a figure  of  about  40  m/sec 
for  the  velocity,  which  is  of  the  order  of  the  velocities  determined  inde- 
pendently by  other  means. 


frequency  200  hZ  before  computing. 
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Figure  6.6.  Separation  of  fixed  and  mobile  probes,  S 
(cm)  vs.  delay  time  for  maximum  cross 
correlation  coefficient,  (milliseconds), 
from  data  of  3/15/77. 


74 


Chapter  7 


SUMMARY  AND  CONCLUSIONS 

This  report  describes  an  experimental  study  of  four-pin  probes  as 
a diagnostic  technique  for  measuring  the  electrical  conductivity  of  flow- 
ing MHD  plasmas.  Emphasis  has  been  placed  on  achieving  temporal  and  spa- 
tial resolution  and  on  investigating  the  accuracy  of  such  probes  in  the 
presence  of  possible  flow  effects. 

After  a brief  account  of  the  theory  of  operation  of  four-pin  probes, 
the  actual  construction  of  the  probes  and  the  relevant  features  of  the 
arcjet  facility  used  to  generate  a flowing  seeded  argon  test  plasma  were 
described.  This  was  followed  by  the  results  of  experimental  investiga- 
tions of  the  i-6v  characteristics  of  the  probes,  in  which  the  linearity 
of  the  characteristics  was  demonstrated  and  the  existence  of  a non-zero, 
fluctuating  inner  pin  voltage  difference,  '^'^offset’  discovered. 

The  presence  of  this  *^^offset  choice,  for  subsequent  experi- 

ments, of  a.c.  probe  excitation  (at  frequencies  in  the  audio  range  well 
above  the  highest  fluctuation  frequencies  present  in  the  flowing  plasma) 
and  the  development  of  the  final  electronic  signal-processing  system.  This 
system  was  finally  employed  in  experiments  to  determine  how  the  indicated 
conductivity  readings  from  four-pin  probes  changed  with  the  flow  velocity 
over  the  pins  and  with  the  orientation  of  the  pins  with  respect  to  the 
flow.  These  experiments  also  included  some  analog  tape  recordings  of  time- 
resolved  conductivity  fluctuations,  from  which  various  statistical  proper- 
ties were  later  deduced. 

The  major  conclusions  which  we  have  drawn  about  four-pin  probes  as 
a result  of  this  research  are  as  follows: 

• The  i-6v  characteristics  of  the  probes  used  remained  linear  if  the 

amplitude  of  the  a.c.  voltage  applied  to  the  outer  pins  was  kept  below 

about  1 volt.  This  criterion  seemed  to  be  a more  useful  rule  of  thumb  for 

ensuring  linearity  than  some  criterion  based  on  the  current  density  at  the 

outer  pin  surfaces,  which  varied  with  probe  geometry  but  was  usually  less 

2 

than  approximately  0.2  amps/cm  . 
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■ Even  when  no  current  passes  between  the  outer  electrodes,  a finite 
and  fluctuating  voltage  difference,  ‘^'^offset’  between  the  inner 

pins.  This  fact  requires  that  the  probe  not  be  operated  with  d.c.  current 
if  erroneous  results  are  to  be  avoided. 

• When  excited  at  frequencies  well  above  the  highest  frequency  of  the 
naturally  occurring  conductivity  fluctuations  in  the  flow,  the  four-pin 
probe  system  described  in  this  report  appears  to  yield  reliable  and  accu- 
rate time-resolved  data  on  conductivity  fluctuations.  That  is,  all  sta- 
tistical parameters  characterizing  the  fluctuations  are  accurate,  except 
possibly  the  mean  value,  O. 

• As  discussed  in  Appendix  A,  the  probe  constant  could  only  be 

determined  to  a precision  of  approximately  5%  to  10%  by  the  method  of  im- 
mersion in  KC2,  solutions.  This  fact  limits  the  probe's  absolute  accuracy 
even  in  the  absence  of  flow  effects. 

• For  a fixed,  "end-on"  probe  orientation,  there  was  a significant  de- 
crease in  the  probe's  indicated  conductivity  as  the  velocity  Increased 
(cf.  Fig.  5.3).  One  way  of  interpreting  these  results  is  in  terms  of  a 
velocity-dependent  correction  factor  to  the  probe  constant  (cf.  Fig.  5.4). 
Based  on  the  present  measurements,  it  seems  that  failure  to  account  for 
this  velocity  effect  could  result  in  an  inferred  conductivity  as  much  as 
two  or  three  times  lower  than  the  actual  conductivity. 

• There  is  also  an  effect  of  the  probe  pin  orientation  relative  to  the 
flow  direction,  which  manifests  Itself,  for  a fixed  flow  velocity,  as  a 
decrease  in  the  indicated  conductivity  as  the  pins  are  rotated  from  the 
"end-on"  configuration  towards  a sidewise  orientation  (cf.  Fig.  5.6). 

• The  following  comments  can  be  made  regarding  the  suitability  of  four- 
pin  probes  as  a practical  diagnostic  for  MHD  work  in  general.  For  rela- 
tive, time-resolved  conductivity  fluctuations  measurements,  such  probes 
would  seem  to  be  very  useful.  (No  experiments  were  performed  to  assess 
the  effect  of  a magnetic  field  on  four-pin  probes,  though  this  would  prob- 
ably complicate  matters  considerably.  Also,  materials  problems  would  have 
to  be  overcome  in  the  case  of  higher-temperature,  oxygen-bearing,  possibly 
slag-laden  combustion  flows.)  For  absolute  conductivity  measurements 
(either  time-resolved  or  averaged) , relatively  large  flow  effects  would 
appear  to  limit  the  accuracy  of  four-pin  probes  from  a practical  standpoint. 
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Using  a particular-sized  and  -oriented  probe,  we  have  found  that  the  flow 
effects  can  be  compensated  for  by  a suitable,  empirically  based  correction 
factor.  To  expect  good  absolute  accuracy  from  any  other  probe,  of  differ- 
ent size,  differently  oriented,  and  in  different  flow  conditions  would  seem 
to  require  that  flow-effect  calibration  experiments  be  carried  out  on  that 
particular  probe  also.  While  this  is  of  course  always  possible,  in  prac- 
tical cases  it  may  not  always  be  feasible  or  worthwhile. 
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Appendix  A 

EXPERIMENTAL  DETERMINATION  OF  PROBE  CONSTANT 


In  this  appendix  we  present  some  detailed  experimental  procedures  which 
were  found  to  be  necessary  in  order  to  obtain  reasonably  accurate  measurements 
of  the  probe  constants,  k^,  for  the  four-pin  probes  used  in  the  research 
reported  here.  The  actual  electronics  setup  used,  plus  the  measurement  and 
data  reduction  procedures  for  determining  have  already  been  presented  in 

Section  4.4.  The  material  in  this  appendix,  on  the  other  hand,  relates  mainly 
to  problems  arising  from  the  use  of  a liquid  electrolyte  (aqueous  solutions  of 
KCl)  as  a calibration  medium. 

It  was  immediately  apparent  from  preliminary  attempts  to  calibrate  probes 
in  KCl  solutions  that  d.c.  excitation  of  the  probe  uas  highly  unsuccessful  due 
to  gas  bubble  formation  on  the  electrodes  and  to  the  large  and  unsteady  local 
concentration  gradients  set  up  in  the  electrolyte  around  the  probe  pins.  Thus, 
while  it  was  clear  that  a.c.  excitation  was  required,  it  was  not  known  at 
what  frequency  the  probe  should  be  excited.  When  k^  was  measured  (as  described 
in  Section  4.4)  over  a range  of  frequencies,  different  results  were  obtained. 

A variation  of  about  30%  was  measured  in  k^  over  the  frequency  range  100  Hz 
to  50  kHz,  for  one  particular  probe  (see  Figure  A-1). 

Further  investigation  into  the  electrochemistry  of  electrolyte  solutions 
revealed  that  similar  problems  had  been  encountered  by  the  original  chemical 
researchers  attempting  to  measure  the  electrical  conductivity  of  electrolytes 
to  high  accuracy.  Most  of  the  observed  frequency  dependence  at  low  frequencies 
was  attributed  to  the  effects  of  polarization  of  the  electrolyte  (such  as 
electrolysis,  causing  gas  bubble  deposits  which  change  the  effective  area  of 
the  electrodes),  while  at  higher  frequencies  (roughly,  above  10  kHz)  errors 
due  to  the  capacitative  effects  of  ionic  sheaths  around  the  electrodes  begin 
to  become  appreciable.  The  standard  frequency  range  for  electrolytic  con- 
ductivity work  was  in  the  range  of  several  hundred  hertz  to  several  kilohertz. 
The  standard  reference  for  reduction  of  polarization  errors  is  a classical 
series  of  papers  in  the  American  Chemical  Society  Journal  (listed  in  reference 
[13])  describing  the  research  done  by  Grinnell  Jones  and  coworkers  at  Harvard 
University  in  the  early  1930's. 
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3/15/77,  before  and  after  platim'zation  of  the  electrodes,  as  a function  of 
excitation  frequency. 


r 


According  to  [l3],  polarization  errors  can  be  greatly  reduced  or 
eliminated  by  making  the  electrodes  inert  and  electrochemically  reversible, 
and  the  standard  technique  for  achieving  this  is  to  electroplate  -the  elec- 
trodes with  very  finely  divided  black  platinum  particles — a process  invented 
by  Kohlrausch  in  1897,  called  platinization.  Following  the  prescription  in 

[13],  all  four  of  the  probe  electrodes  were  platinized  by  passing  a total 

2 

charge  of  at  least  1.7  coulombs  per  cm  of  electrode  area  (at  current  levels 
about  10  ma)  through  each  electrode  while  immersed  in  a platinizing  solution 
of  0.025  Normal  hydrochloric  acid,  0.3%  by  weight  platinum  chloride,  and 

0.025%  by  weight  lead  acetate  (The  actual  solution  contained  1 gram  PtCl  , 

3 3 ^ 

333  cm  H^O,  0.255  cm  HCl,  and  0.083  gm  lead  acetate).  After  platinization, 

the  black-coated  pins  were  very  carefully  washed  in  distilled  water,  ana  then 
immersed  in  the  standard  1 Normal  KCl  solution  (prepared  by  diluting  74.555  gm 
KCl  weighed  in  air  to  1 liter  total  volume  with  distilled  water);  the  con- 
ductivity of  this  electrolyte  is  given  as  a function  of  temperature  in 
reference  [14]. 

In  addition  to  platinization  of  the  electrodes,  several  other  more  routine 

precautions  were  taken.  For  each  measurement  of  k^ , a beaker  of  "fresh" 

KCl  solution,  taken  from  the  same  carefully  prepared,  sealed  reserve  supply, 

was  used,  and  the  beaker  was  sealed  off  to  avoid  evaporation  of  H^O  after  the 

probe  was  immersed  in  the  solution.  An  accurate  mercury-in-glass  thermometer 

with  0.1°C  graduations  was  used  to  monitor  the  solution  temperature.  While 

making  the  electrical  measurements,  current  was  allowed  to  flow  only  long 

enough  (typically  about  2 seconds)  to  achieve  (and  read  on  a digital  voltmeter) 

a steady  divider  output  voltage,  and  for  all  measurements  the  current  density 

2 

was  held  to  the  same  value  (=  0.2  amps/cm  at  pin  surfaces)  by  changing  the 
amplitude  of  the  applied  a.c.  excitation  voltage  as  required. 

With  platinization  (and  the  above  precautions  heeded),  the  measured  values 
of  k^  showed  a less  severe  frequency  dependence  than  was  obtained  without 
platinizing  the  electrodes.  A comparison  of  the  measured  frequency  dependence 
of  k for  one  particular  probe,  before  and  after  platinization,  is  shown  in 
Figure  A-1.  The  variation  of  k witli  frequency  was  never  successfully 
eliminated  entirely  on  any  of  the  probes  used,  and  there  was  some  variability 
of  the  severity  of  the  frequency  dependence  from  probe  to  probe.  Generally 
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the  value  of  k at  1 kHz  was  used  to  reduce  the  experimental  data  on  plasma 

P . ... 

conductivity,  but  the  estimated  uncertainty  in  this  value  could  be  as  much  as 
5%  to  10%,  depending  on  the  probe.  The  curves  of  k^  versus  frequency  for  the 
particular  fixed  and  mobile  probes  used  in  the  most  successful  experiments 
(upon  which  the  important  results  of  Figure  5.3  are  based)  are  shown  below  in 
Figure  A-2.  There  was  about  a 6%  reduction  relative  to  the  pre-test  value  in 
kp  of  the  mobile  probe  after  it  had  been  used  in  the  experiment,  perhaps 
reflecting  some  small  dimensional  changes. 
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Appendix  B 

REDUCTION  OF  CONDUCTIVITY  FLUCTUATIONS 
INTRODUCED  BY  ARCJET  POWER  SOURCE 


In  the  early  experiments  using  four-pin  probes  in  the  a’^cjet  facility, 
oscilloscope  photographs  of  the  inner  pin  voltage  difference  vs.  time  showed 
a repeatable,  characteristic  pattern  of  very  sharp  peaks  at  a frequency  close 
to  360  Hz.  These  fluctuations  were  interpreted  as  indicating  the  presence  of 
sharply  peaked  electrical  conductivity  fluctuations.  It  was  hypothesized 
that  these  o fluctuations  were  due  to  360  Hz  temperature  fluctuations 
probably  caused  by  the  "ripple"  component  of  the  current  supplied  to  the  arcjet 
by  the  laboratory's  d.c.  power  supply.  This  appendix  describes  how  these 
fluctuations  were  greatly  reduced  by  employing  the  large  water-cooled  MHD 
magnet  as  a inductor  in  series  with  the  arcjet. 

The  d.c.  power  source  consists  of  twelve  arc  welder  power  supplies  (A.O. 
Smith  Co.,  Model  10,000A)  connected  as  shown  in  Figure  B-la.  The  basic  wiring 
diagram  for  each  of  these  individual  units  is  given  in  Figure  B-lb.  This 
system  produces  a three-phase,  full-wave  rectified  output  voltage  waveform 
whose  shape  is  approximately  that  obtained  by  linearly  superimposing  three 
full-wave  rectified  sinusoidal  waveforms  shifted  in  phase  from  each  other  by 
120°,  which  yields  a d.c.  output  with  a 360  Hz  "ripple"  component  whose 
amplitude  is  several  percent  of  the  d.c.  value. 

The  original  connection  of  this  pov/er  supply  to  the  arcjet  is  diagramed 
in  Figure  B-2.  As  can  be  seen  by  examination  of  Figure  B-2  and  the  circuit 
diagrams  of  Figure  B-1,  the  current  which  passes  through  the  arcjet  negotiates 
the  complete  circuit  without  encountering  resistive  or  inductive  curcuit 
elements  which  would  tend  to  impede  current  fluctuations.  Furthermore,  the  V-I 
characteristics  of  a heavy,  high-pressure  d.c.  arc  such  as  exists  in  the  arcjet 
are  almost  completely  flat  (cf.  Reference  [12],  p.365);  that  is,  the  voltage 
drop  is  determined  by  the  gap  length,  the  gas  composition  and  the  gas  flow  rate 
through  the  arc,  and  is  essentially  independent  of  the  current  level,  once  the 
arc  has  been  struck.  This  fact  means  that  the  arc  itself  will  not  offer  any 
resistance  to  fluctuations  in  the  current.  Consequently  the  current  through 
the  arcjet  will  follow  the  fluctuating  eraf  of  the  power  supply  units  and  will 
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have  a fluctuating  component  at  360  Hz  which  is  several  percent  of  the  mean 

current;  this  current  fluctuation  would  produce  a small  360  Hz  gas  temperature 

fluctuation,  which  would  lead  to  a significant  360  Hz  conductivity  fluctuation. 

The  obvious  solution  to  this  problem  was  to  place  a large  inductor,  one 

3 

capable  of  carrying  the  full  arc jet  current  (of  the  order  of  10  amperes),  in 
series  with  the  arc jet.  The  inductive  impedance  u)L  would  then  hopefully 
be  large  enough  to  satisfactorily  attenuate  the  current  fluctuations. 
Fortunately,  such  an  inductor  was  readily  available  in  the  laboratory  in  the 
form  of  a large  2.7  tesla  iron-core,  water-cooled  electromagnet  used  for  MHD 
experiments . 

The  inductance  L of  this  magnet  was  measured  using  the  simple  setup 
shown  in  Figure  B-3.  The  time  response  V(t)  of  the  voltage  across  the 
known  resistor  R was  recorded  on  an  oscilloscope  after  closing  switch  S. 
Since 

V(t)  = Vo^  1 - e ^ ^ ^ 

a plot  of 


vs.  t will  be  a straight  line  with  slope  equal  to  R/L.  The  value  of  L 
obtained  in  this  fashion  was  L = 0.74  henry  with  the  pole  pieces  fully  closed. 
The  impedance  wL  offered  by  this  magnet  to  360  Hz  fluctuations  is  thus 
about  1675  ohms,  whereas  the  impedance  offered  to  the  steady  d.c.  current 
is  only  0.214  ohms,  which  is  the  d.c.  resistance  of  the  magnet  coils. 

The  magnet  (as  an  inductor)  was  connected  in  series  with  the  arcjet  as 
shown  in  Figure  B-4.  The  diode  shown  in  this  figure  is  a Motorola  type 
MR1269FL  silicon  power  rectifier  capable  of  carrying  the  full  arcjet  current 
in  the  forward  direction.  In  the  event  of  a sudden  extinguishing  of  the  arc 
for  whatever  reason,  this  diode  protects  the  magnet  coils  from  internal  arcing 
due  to  high  L di/dt  voltages  by  allowing  the  stored  current  to  smoothly  decay 
to  zero;  the  stored  energy  is  then  dissipated  in  the  ohmic  resistance  of  the 
water-cooled  magnet  coils. 
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Figure  B-4.  Connection  of  MHD  magnet  in  series  with  the  arcjet  as  an  inductor 
to  reduce  current  fluctuations  introduced  by  the  laboratory  d.c. 
supply  system. 
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In  addition  to  the  large  diode,  there  are  five  small  "Transzorbs" 

(trade  name  for  General  Seminconductor  Industries  type  KEIOOCA  bipolar  transient 
suppressors)  connected  in  series  with  one  another  and  then  in  parallel  across 
the  diode  (cf.  Figure  B-4).  These  effectively  short  out  any  high-speed,  low- 
energy  transient  voltage  pulses  in  excess  of  500  volts,  and  were  installed  as 
a precaution  to  protect  both  the  large  diode  and  the  magnet  windings  from 
possible  damage  from  the  high  voltage  arc  starter. 

The  system  as  shown  in  Figure  B-4  performed  entirely  satisfactorily;  no 
malfunctions  occurred,  and  the  arcjet  current  ripple  fraction  was  reduced  to  a 
measured  value  of  0.01%. 
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DETERMINATION  OF  THE  FLOW  VELOCITY 
ALONG  THE  NOZZLE  AXIS 
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The  nozzle  employed  in  these  experiments  was  machined  from  a solid 
cylinder  of  boron  nitride,  and  the  actual  contour  was  chosen  somewhat 
arbitrarily  at  the  time  of  machining  to  provide  a smooth,  slowly  converging 
shape.  Figure  C-1  is  an  approximately  correctly  proportioned  sketch  showing 
the  geometrical  details  of  the  flow  path  in  the  vicinity  of  the  nozzle.  The 
maximum  mean  velocitv,  U , occurs  at  the  minimum-area  section  at  the  exit, 
and  is  given  approximately  by  the  formula 


U 

max 


R m T 


(C-1) 


where  R is  the  universal  gas  constant,  M the  molecular  weight  of  the  gas, 

m the  mass  flow  rate,  P the  pressure,  d . the  minimum  diameter  of  the 

min 

nozzle,  and  T the  gas  temperature.  For  argon  at  one  atmosphere  pressure 

and  2200°K  (a  typical  estimated  gas  temperature  in  these  experiments),  and 

with  m = 0.03  Ibm/sec  and  d . = 1.1  inch  (used  in  all  experiments),  formula 

min 

(C-1)  yields  a value  U = 100  m/sec.  This  value  leads  to  a Mach  number 

max 

less  than  about  0.1,  which  means  that  for  purposes  of  determining  the  flow 

velocities  along  the  nozzle  axis,  negligible  error  is  made  by  assuming  the 

flow  to  be  incompressible.  The  remainder  of  this  appendix  describes  briefly 

the  three  methods  used  to  estimate  U(z)/U  as  a function  of  the  distance  z 

max 

along  the  axis  from  the  exit  plane.  The  resulting  three  curves  of  U(z)/U 

max 

vs.  z are  shown  superimposed  in  Figure  C-2. 

The  first  and  simplest  method  was  to  make  the  approximation  that  the 
flow  (in  addition  to  being  inviscid  and  incompressible)  was  strictly  one- 
dimensional, in  which  case  we  obtain  from  the  continuity  condition  the  result 


U(z)/U  - f d . /d(z)V  . (C-2) 

max  y min  / 

The  diameter  d(z)  was  carefully  measured  as  a function  of  z by  mounting 
the  nozzle  on  the  table  of  a vertical  milling  machine  and  tracing  the  contour 
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Figure  C-2.  Centerline  velocity  ratio  u/umax  ''s.  axial 

distance,  z,  into  nozzle  from  the  exit  plane, 
as  determined  by  three  methods. 
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with  a dial  indicator  while  recording  the  vertical  and  horizontal  translations 
of  the  table  required  to  keep  the  indicator  in  contact  with  the  inner  surface 
of  the  nozzle.  This  one-dimensional  calculation  will  obviously  give  erroneous 
results  whenever  there  is  a sudden  discontinuity  in  the  area,  as  occurs,  for 
example,  at  the  nozzle  entrance  (see  Figure  C-2). 

The  second  method  was  to  solve  for  the  two-dimensional  velocity  field 
on  a computer,  thereby  improving  the  accuracy  of  the  solution,  with  respect 
to  the  previous  one-dimensional  method.  The  equation  solved  was 


_ i li  + Ijii  = n 

.2  r 3r  .2  " ' 

3r  3z 


(C-3) 


where  is  the  Stoke 's  stream  function.  This  equation  governs  the 

axisymmetric , irrotational  flow  of  an  inviscid,  incompressible  fluid.  The 
radial  and  axial  velocity  components  are  found  from  the  solution  according 


to  the  formulas 


r 3z  * 


= i 

r 3r 


(C-4) 


The  boundary  condition  employed  was  to  specify  i//  all  around  the  domain 
boundary;  this  domain,  chosen  to  approximate  resonably  closely  the  actual  flow 
geometry,  was  as  shown  in  Figure  C-3.  Also  shown  are  the  plotted  streamlines 
obtained  by  the  computer  code  "AXI.AX",  which  was  generously  provided  by 
Mr.  Marion  Jenkins  of  the  High  Temperature  Gasdynamics  Laboratory  at  Stanford. 
The  stream  function  ij;  was  chosen  to  be  constant  along  the  centerline  and 
along  the  channel  walls.  The  velocity  profiles  at  the  inlet  and  outlet  were 
assumed  to  be  uniform  (v  constant,  v = 0),  so  that  the  values  of  specified 

Z IT  2 

along  the  inlet  and  outlet  had  to  satisfy  the  condition  '?  = (constant)  r , 
obtained  by  integrating  the  second  of  equations  (C-A). 

As  seen  from  Figure  C-2,  the  results  of  the  two-dimensional  numerical 
solution  are  in  fairly  close  agreement  with  the  one-dimensional  results  except 
near  the  ends  of  the  nozzle,  where  the  former  results  appear  more  realistic. 

The  main  limitations  of  the  numerical  method  are  its  failure  to  account  for 
visous  effects,  compressibility  effects,  the  presence  of  swirl  (non-zero  V0), 
the  possibly  non-uniform  inlet  and  exit  velocity  profiles,  and  the  lack  of 
complete  geometrical  identity  between  the  actual  flow  region  and  the  domain 
used  for  the  numerical  solution. 


Due  to  uncertainty  regarding  the  quantitative  effects  of  the  limitations 

just  mentioned  on  the  calculated  axial  velocity  ratio  U(z)/U  , it  was 

^ max 

decided  to  obtain  a direct  experimental  determination  of  this  ratio  using  a 
Pitot  pressure  probe.  The  local  velocity  is  found  from  the  formula 


where  p and  p are  the  stagnation  pressure  and  density,  respectively,  and 
o o 

P is  the  local  static  pressure.  Since  the  flow  can  be  closely  approximated 

as  incompressible,  we  may  divide  by  the  corresponding  formula  for 

obtain 


U(z) 

U 

max 


r(p  -p(z) 

o 

p -p  . ) 
o min 


(C-5) 


The  pressure  differences  were  measured  with  a sensitive  inclined  manometer  set 
up  as  shown  in  Figure  C-A.  The  pressure  P(z)  was  sensed  by  a probe  con- 
structed from  3/32-inch  diameter  tubing  hemispherically  capped  and  drilled 
with  four  0.020  inch  holes  at  a section  1.25  inches  downstream  from  the  capped 
end.  Hiis  probe  was  mechanically  traversed  along  the  axis  of  the  nozzle 
while  room  temperature  nitrogen  was  flowing  through  the  plenum  and  nozzle  with 
a maximum  Mach  number  of  0.1. 

The  experimental  results  for  ^ obtained  from  this  static 

pressure  measurement  technique  are  shown  on  Figure  C-2,  and  indicate  a slightly 

more  gradually-changing  velocity  near  the  nozzle  entrance  and  exit  than  is 

predicted  by  the  previously-described  calculational  methods.  Tlie  error  bars 

shown  in  Figure  C-2  reflect  mainly  the  uncertainty  in  reading  AP  from  the 

manometer,  whose  fluid  level  was  fluctuating  somewhat  during  the  measurements. 

Due  to  the  relative  "directness"  of  this  experimental  technique  for  determining 

U(z)/U  , as  opposed  to  the  calculational  techniques  which  rely  on  many 
max 

assumptions,  the  experimentally-determined  curve  of  U(z)/U  vs.  z was 

max 

deemed  more  accurate  and  was  used  in  all  subsequent  data  presentations. 
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Appendix  D 

DETAILS  OF  FOUR-PIN  PROBE 
EXCITATION  AND  SIGNAL  PROCESSING  CIRCUITRY 


This  appendix  consists  of  a set  of  wiring  and  circuit  diagrams  for  the 
probe  excitation  circuitry  and  for  the  signal  conditioning  equipment  employed. 

The  probe  excitation  circuitry  consisted  of  two  units:  a floating 
function  generator (constructed  from  a kit,  shown  in  Figure  D-2)  to  provide  a 
low-power  excitation  voltage,  and  a floating  power  amplifier  (shown  in  Figure 
D-1)  designed  to  amplify  the  excitation  signal  supplied  by  the  function 
generator  to  power  levels  sufficient  to  drive  reasonably  large  currents 
between  the  outer  probe  pins.  This  combination  of  equipment  permitted  exciting 
the  probes  over  the  frequency  range  0.1  Hz  to  around  100  kHz. 

A single  unit  (the  "signal  conditioning  unit")  was  constructed  to 
accomplish  most  of  the  signal  processing  vi/hich  is  indicated  in  block  form  in 
Figure  4.4  and  4.6.  After  this-  signal  conditioning  unit  was  constructed, 
several  additional  commercial  differential  amplifiers  and  electronic  filters 
were  acquired,  and  these  were  used,  together  with  several  of  the 
circuits  in  the  existing  signal  conditioning  unit,  to  form  the  final  elec- 
tronics setup  used  in  the  experiments.  This  final  setup  is  shown  schematically 
in  Figure  D-3. 

Figure  D-4  shows  the  chassis  wiring  of  the  signal  conditioning  unit 
itself.  The  remaining  figures,  D-5  through  D-10,  are  detailed  circuit 
diagrams  of  the  individual  units  (differential  amplifier,  low-pass  filter, 
precision  linear  rectifier,  fixed-  and  variable-gain  amplifier,  and  polarity 
inverter)  which  are  contained  within  the  signal  conditioning  unit.  The 
analog  divider  is  mounted  in  this  unit  also,  but  since  it  is  a commercially 
available  device  (Analog  Devices,  Inc.,  Model  436B,  Two-Quadrant  Divider), 
no  circuit  diagram  is  included  here. 


LH340-K 


Dur-pin  probes. 


Figure  D-3.  Schematic  drawing  of  final  hookup  of  Signal  Conditioning 
Unit  and  external  commercial  differential  amplifiers  and 
filters,  to  form  the  electronics  configuration  of 
Figure  4.4. 


Figure  D-5.  High  Input  Impedance  Unity-Gain  Differential  Amplifier 
Circuit.  Nominal  specifications:  bandwidth  = Mhz, 
CMRR  =.  50,000  at  10  Hz,  common  mode  range  = + 11  volts 
input  impedance  = 10 ohms,  output  impedance  = 1 ohm, 
output  swing  = ± 10  volts. 
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Figure  D-6.  Low-Pass  Filter  Circuit.  Filter  type  is  multiple  feedback 
fourth-order,  with  Butterworth  (-80  db/decade  rolloff) 

response,  f^  = -3  db  cutoff  frequency.  Gain  = 1 + R5/R5 
(=2.0  for  values  shown). 
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Appendix  E 

PRESERVATION  OF  LOU-FREQUENCY 
CONDUCTIVITY  FLUCTUATIONS 


The  argument  to  be  presented  here  is  intended  to  demonstrate  that  the 

variations  in  slope  of  i-6V  characteristics  such  as  these  presented  in 

Figure  4.2,  obtained  as  explained  in  Section  4.1  with  the  high-pass  filter 

setup  of  Figure  4.1,  are  in  fact  directly  proportional  to  the  electrical 

conductivity  fluctuations,  even  though  the  latter  nay  occur  at  frequencies 

far  below  the  cutoff  frequencies  of  the  filters. 

To  begin  with,  we  return  to  the  analysis  of  Section  1.3  and  write  down 

expressions  for  the  instantaneous  current  i(t)  and  inner  pin  voltage 

difference  6V.  (t).  Referring  to  the  notation  in  Figure  1.2  and  applying 

inner  “ ” 

Kirchoff's  law  around  the  complete  circuit,  we  obtain 

V . - i R - AV  + AV  T 

applied  sh  WF,3  sheath, 3 31 

-k-,pi  - AV  , .,  , + AV,_  , = 0 , 

24  sheath, 4 WF,4 


which,  upon  rearrangement,  may  be  written 


V - . .(t)  - A V , j (t) 

appl  led anode-cathode 

i(t)  = ^ ^ot 


(E-1) 


In  equation  (E-1)  p(t)  is  the  plasma's  resistivity,  and  ^'^gnode-cathode 

k are  defined  as  follows: 
tot 

4V  . .u  j - - AV^  ,)  - (AV  . - AV  , (F-2) 

anode-cathode  VfF,3  WF,4  sheath, 3 sheath, 4 


k.  ^ = k.,  + k + k . . 
tot  31  p 24 


(E-3) 


Here  k is  constant  in  time  (since  it  is  a function  of  the  geometry  of 
tot 

the  probe),  but  ^'^anode-cathode  ® time-varying  quantity  analagous  to 
^^offset’  addition  to  (E-1)  we  have  directly  from  equation  (1.3)  that 

6V.  = k p i + iSV  . , 

inner  p offset 
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which  becomes,  after  substitition  of  the  result  (E-1)  for  i, 


6V.  (t)  = 6V  --  ^(t)  + k p( 

inner  offset  p 


«>[ 


^applied^*^^ 


AV 


anode-cathode 


R . + k p (t) 

sh  tot 


(t)] 


(E-4) 


Now  we  consider  the  outcome  of  passing  the  i and  6V.  signals 

given  by  (E-1)  and  (E-A)  through  the  high-pass  filters  sho^m  in  Figure  4.1. 

The  transfer  function  of  these  filters  is  represented  schematically  in 

Figure  E-1,  in  which  the  relative  magnitudes  of  the  relevant  frequencies  is 

indicated  on  the  frequency  axis.  From  this  figure  we  see  that  only  the 

"fast-varying"  signals  whose  frequencies  are  close  to  f^  (typically  20  kHz) 

or  greater  will  be  transmitted  by  the  filter.  Keeping  in  mind  that  the 

quantities  AV  , , , (t),  6V  ..(t)  and  p(t)  are  "slowly-varying," 

anode-cathode  ’ offset 

at  least  relative  to  V .(t),  we  can  assert  (and  could  demonstrate  more 

applied 

formally  using  Fourier  series)  that  the  output  signals  will  contain  only  those 

terms  in  which  V ,(t)  is  a factor,  i.e., 

applied  ’ 


^filtered 


V , . .(t) 

applied 

R , + k p(t) 

sh  tot 


(E-5) 


and 


6V.  (t)  = k p(t) 

inner  p 

filtered 


r V , . 

appl  1 

R + k 
sh 


ed 

k^  ^p(t) 
tot 


(E-6) 


Substituting  (E-5)  into  (E-6),  and  replacing  the  resistivity  p(t)  by  the 
inverse  of  the  conductivity  o(t)  we  obtain 


i,._  ,(t)  = 6v.  ...  .(t)  , 

filtered  k inner, filtered 


(E-7) 


which  shows,  as  asserted  above,  that  the  fluctuating  slope  of  the  i-6V 
characteristics  (obtained  after  high-pass  filtering  as  in  Figure  4.1)  is 
proportional  to  o(t). 
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Figure  E-1.  Approximate  transfer  function  of  high  pass  filters  used  in 
4-pin  probe  signal  processing. 
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Appendix  F 

ESTIMATION  OF  POSSIBLE  DECREASES  IN 
CONDUCTIVITY  ALONG  THE  NOZZLE 


One  of  the  major  experimental  results  of  the  present  research  was  a 
measured  decrease  of  around  50%  in  the  electrical  conductivity  indicated  by 
a mobile  4-pin  probe  as  the  probe  was  moved  along  the  axis  of  the  subsonic 
nozzle  from  the  upstream,  low-velocity  region  (Station  1)  to  the  downstream, 
high-velocity  region  (Station  2).  We  discuss  in  this  appendix  to  what 
extent  the  "true"  local  conductivity  may  have  changed  between  these  two 
stations  due  to  naturally  occurring  processes  in  the  plasma. 

To  begin  the  discussion,  we  review  the  known  facts  about  the  state  of 
the  gas  at  station  1 (in  the  plenum  region).  During  the  experiment  which 
yielded  the  results  of  Figure  5.3,  the  tungsten-rhenium  thermocouple  in  the 
plenum  core  indicated  a junction  temperature  of  about  1700  K.  Using  various 
estimates  of  the  radiation  correction  to  this  thermocouple  reading,  we 
arrive  at  the  result  that  the  heavy  gas  temperature  Tg  lies  approximately 
in  the  range  of  1900  K to  2100  K.  The  electrical  conductivity  indicated  by 
both  the  fixed  plenum  probe  and  the  mobile  probe  at  the  low-velocity  station  1 
was  13  ± 1 mho/m.  From  calculations  of  electron  transport  properties  by 

Brown  (specifically.  Appendix  2.4  of  Reference  6)  we  can  then  infer  that 

12-3  . — 

N = 4.3  X 10  cm  and  the  collision  frequency  v is  approximately 

® 9 -1  — . 

9.3  X 10  sec  (v  „ is  insensitive  to  T over  the  range  1400  K to  2200  K). 
eH  e ® 

We  now  wish  to  know  by  how  much  a (or  N^,  which  is  directly  propor- 
tional to  o in  our  case)  might  change  as  the  gas  flows  down  the  nozzle.  In 
general,  this  is  a difficult  problem.  To  the  extent  that  the  electron  gas 
is  only  weakly  coupled  to  the  heavy  particle  gas,  one  could  in  principle 
answer  this  question  by  solving  simultaneously  suitable  forms  of  the  electron 
energy  equation  and  the  electron  continuity  equation,  which  constitute  a 
pair  of  coupled,  first-order  differential  equations  for  N^  and  T^  as 
functions  of  time.  We  do  not  propose  such  a calculation  here.  Rather,  we 
ask  what  would  be  the  possible  limiting  cases  if  such  detailed  calculations 
were  carried  out. 
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One  possible  limiting  case  (which  we  label  "Case  A")  would  be  that  in 
which  local  thermodynamic  equilibrium  (i.e.  T = T ) is  maintained  through- 
out the  flow.  In  this  case  we  would  need  to  examine  mechanism  leading  to 
decreases  in  T , and  then  relate  the  magnitudes  of  such  decreases  AT  to 
changes 


AN 


Ao 

o 


in  N and  o via  the  Saha  equilibrium  relationship  between  N and  T . 
e e e 

Another  possibility  (which  we  call  "Case  B")  would  be  that  Saha 
equilibrium  is  maintained  (in  the  sense  that  N and  would  be  related 

through  the  Saha  equation)  at  each  point  in  the  flow,  but  that  T T (in 
fact  T would  be  less  than  T for  our  experimental  conditions).  In  this 

e g _ , 

case  it  would  be  necessary  to  examine  mechanisms  leading  to  decreases  in  T , 

e 

such  as  radiation  escape  from  the  plasma,  and  to  compute  the  corresponding 
changes  in  o from  the  Saha  equation.  It  is  easy  to  show  that  the  relation- 


ship between  the  relative  fractional  changes  in 
these  two  cases  is 


do  _ e_  _ + 1 1 

o “ Ne  U 2kT^ 


T 


and  T in  either  of 


(F-1) 


where  is  the  seed  ionization  potential  and  T is  either  T or  T , 

. . eg’ 

whichever  is  applicable. 

At  the  other  extreme  ("Case  C")  lies  the  possibility  that  the  electron 

density  (and  hence  the  conductivity)  is  essentially  frozen  at  its  initial 

value.  This  condition  could  occur  if  the  electron-ion  recombination  time 

(t  , ) were  much  greater  than  the  flow  time  of  the  plasma  down  the  nozzle, 
recomb 


or  T which  might  be  simultaneously 
e 


In  this  case,  any  changes  in  T 
taking  place  would  be  irrelevant  as  far  as  estimating  changes  in  o is 
concerned,  since  o would  be  essentially  constant  dovm  the  nozzle. 

We  consider  first  this  latter  possibility.  The  time  required  for  a 
fluid  element  to  flow  from  station  1 to  station  2 in  these  experiments  was 
determined  from  the  relation 


flow  U 


max 


, \ max/ 
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Using  Figure  C-2  to  provide  the  measured  profile  of  ® function 


of  z,  and  carrying  out  the  integration  graphically,  we  obtained  t 
-3 

1.37  X 10  sec. 

To  estimate  x 

electron-ion  recombination  as 


flow 


, we  write  the  rate  of  change  of  N due  to 
recomb.  e 


2 

TT  = -“^e 


(F-2) 


and  for  the  recombination  coefficient  a we  may  employ  the  Hinnov  and 
Hirschberg  formula  (see  Reference  [8],  Chapter  IX,  Section  4) 


T -9/2 

a = 3.4  X 10  (l0§0^  cm^'sec. 


(F-3) 


Here  is  expressed  in  cm  ^ and  in  °K.  Substituting  (F-3)  into 

(F-2)  leads  to  a differential  equation  which  may  be  integrated  between  t=0 
and  t=T  if  we  assume  to  remain  constant  or  change  slowly.  The  result 


^ [n  ^(t)  N ^(0)] 


= 3.4  x 10 


-99  T -9/2 

nooo-'  * 


(F-4) 


e e 

If  we  somewhat  arbitrarily  define  a characteristic  recombination  time 

T L as  the  time  necessary  for  N to  decrease  to  one  half  of  its  initial 
recomb  e 

value,  we  obtain  from  (F-4)  that 


T 9/2 

f— 

^1000-* 


recomb  2 ^ 10"^^1  N "(0) 

J ® 


(F-5) 


12  -3 

The  initial  was  (as  previously  noted)  4.3  x 10  cm  . There  was  no 

direct  measure  of  T , so  we  use  the  value  T = 1650  K,  arrived  at  as 
explained  below  in  a calculation  taking  radiation  escape  into  account. 
Substituting  these  values  of  and  T^  into  (F-5)  yields  the  result 

-3 

T , = 2.27  X 10  sec. 
recomb 
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Comparing  the  two  characteristic  times  we  find  that  x = 

recomb  flow 

1.66,  or  roughly  that  x u ^ . If  it  had  turned  out  that  x . » 

° recomb  flow  recomb 

X,,  , one  would  be  assured  of  frozen  flow  (Case  C),  while  if  we  had  found 

flow’ 

X , <<  x^,  , we  could  confidently  assume  equilibrium  conditions  (Case  A 

recomb  flow 

or  Case  B)  to  hold.  What  seems  likely,  therefore,  is  that  (and  o)  is 
partially  frozen — decreasing  somewhat  as  the  gas  travels  down  the  nozzle,  but 
less  than  it  v/ould  decrease  if  equilibrium  were  maintained.  The  estimates 
of  conductivity  changes  which  are  made  below  assuming  either  of  the  equi- 
librium situations.  Case  A and  Case  B,  may  thus  be  regarded  as  upper  bounds 
on  the  actual  decrease  in  o. 

Before  proceeding  to  such  calculations,  we  remark  that  the  notion  that 

the  conductivity  may  have  been  frozen,  at  least  partially,  appears  to  be 

supported  by  the  experimental  evidence  of  Cool  and  Zukoski  [ 15  ].  Under 

similar  experimental  conditions  (relatively  slowly  moving,  potassium-seeded 

argon  heated  to  2000  K in  an  arcjet)  they  observed  (see  especially  their 

Figure  4)  that  for  low  applied  current  densities  the  conductivity  leveled 

off  at  a value  approximately  corresponding  to  the  equilibrium  conductivity 

2 

at  the  gas  temperature.  They  also  concluded  that  for  J ~ 0.3  amp/cm  the 
approximate  "steady  state"  form  of  the  electron  energy  equation. 


2 2 • • 

oE  = J /o  = 0 + R , 


(F-6) 


combined  with  the  Saha  equation  at  the  electron  temperature,  offers  an 
accurate  way  of  predicting  the  electrical  conductivity.  For  low  values 
of  J,  however,  (F-6)  appears  (when  compared  to  the  data  of  Cool  and 
Zukoski)  to  overestimate  the  effect  of  radiation  in  the  electron  energy 
balance.  This  result  would  be  in  accord  with  the  notion  that  a frozen 
electron  density  was  being  observed  as  J 0,  in  which  event  the  use  of 
the  Saha  equation  would  not  be  justified. 

We  now  proceed  to  Case  B.  In  the  present  experiments  J = 0,  so  that 
the  approximate  energy  balance  equation 


zm  _ 

-U  ‘ I — - V N2kfT-T\=R 
i:nL  eh  62  \g  e/ 


(F-7) 
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may  be  used  to  estimate  the  change  in  T between  the  two  stations  as  a 

e 

result  of  radiation  escape.  The  calculation  procedure  is  summarized  in  what 

follows.  Since  N is  uniquely  determined  as  a function  of  T • by  the 
e e 

Saha  equation  (which  we  are  now  assuming  to  be  applicable),  and  since  T 

changes  only  very  slightly  down  the  nozzle  (a  fact  which  will  be  shown  to  be 

true  later),  we  may  treat  T as  a constant  and  conclude  that  the  left-hand 

8 

side  of  (F-7)  is  a function  only  of  T , and  may  be  displayed  graphically  as 
such.  The  right-hand  side,  R,  which  represents  the  radiated  energy  lost 
from  the  gas  per  second  per  unit  volume,  may  be  estimated  using  the  results 
of  McGregor  [111  for  resonance  radiation  escape  from  potassium-seeded  argon. 
From  page  148  of  [ 11], 

R = R'  + n^  R"  , (F-8) 

where  n^  is  the  ground  level  number  density  of  potassium  atoms  and  R' 

and  R"  are  coefficients  calculated  by  McGregor  and  tabulated  in  Table  A. 2 

of  [ 11  ].  In  order  to  "enter"  this  table,  the  input  parameters  T , N , and 

e e 

6 must  be  specified.  Since  we  are  assuming  Saha  equilibrium  holds,  N is 

e 

determined  once  T is  specified.  Here  B is  the  radiation  escape  para- 
e 

meter  whose  definition  and  calculation  are  presented  in  Appendix  H of  [11], 
which  gives  the  formula 


B(4P-4S)  = (1.374  X 


lO"^) 


I / 10~^\  ' I 

\ 1000/  \ yr(cm)  J 


(F-9) 


Here  r is  the  radius  of  the  radiating  plasma,  which  is  assumed  to  be  a 
uniform,  infinitely  long  cylinder. 

The  procedure  for  determining  T^  at,  say,  station  1 is  thus  as 
follows:  First  one  picks  a gas  temperature  T and  plots  a graph  of  the 
left-hand  side  (-51)  of  (F-7)  as  a function  of  T^.  Then,  using  the  nozzle 


radius 
meter  6 


r^  at  station  1 in  formula  (F-9),  one  computes  the  escape  para- 


^ at  station  1.  Next,  using  this  B^ 


and  some  assumed  value  for 


T (and  the  corresponding  "Saha"  value  of  N^);  one  enters  McGregor's  Table 
A. 2 (liberally  interpolating  as  necessary)  and  finds  R'  and  R",  and 
finally  computes  R from  (F-8).  This  process  is  repeated  for  several  more 
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assumed  values  of  T until  enough  points  are  obtained  to  plot  R = R(t  . r,) 

e , e i 

as  a function  of  T on  the  same  graph  on  which  -n  is  plotted.  The 
e 

abscissa  of  the  point  of  intersection  of  the  two  resulting  curves  then  is  the 

value  of  X which  satisfies  the  energy  balance  relation  (F-7)  at  station  1. 
e 

When  this  process  was  carried  out  using  ••  2000  K and  the  remaining  para- 
meters appropriate  to  our  plasma  conditions,  the  value  T = 1646  ± 2 K 

e , 1 

was  obtained.  Then  the  entire  process  was  repeated,  changing  only  the  radius 

r to  be  that  at  station  2,  with  the  result  that  T ^ = 1632  ± 2 K. 

e,2 

Substituting  4T  =T  , -T  „ =14K  for  dT  into  (F-1),  and  using  the 

e e,l  e,2  e 

value  T = 1640  K,  we  obtain  a maximum  percentage  reduction  in  a of 
e 

do/o  = 14%,  considerably  less  than  the  observed  50%  drop  in  indicated 
conductivity. 

Finally  we  consider  Case  A,  in  which  it  is  to  be  assumed  that  T = T 

e g 

all  along  the  nozzle.  One  mechanism  leading  to  a decrease  in  T is  the 

8 ^ 

gasdynamic  cooling  which  occurs  as  the  thermal  energy  of  the  gas  is  shifted 
to  directed  kinetic  energy  as  the  flow  accelerates.  The  estimate  of  the 
magnitude  of  this  effect  is  based  on  the  well-known  formula  (see,  for 
example.  Reference  10)  for  the  radio  of  the  stagnation  temperature,  T^,  to 
the  static  temperature  T in  the  adiabatic  flow  of  a perfect  gas 

T 

2-=  1 ^2 


= 1 + ^ 


where 


(F-10) 

and  M is  the  Mach  number.  Using  k = 5/3  for  argon. 


k = c /c 

p V 

plus  the  maximum  Mach  number  M=0.095  obtained  in  these  experiments,  we 
obtain  from  (F-10)  the  result  that 

T -T 

_2 ^ = 0.003, 

o o 

which  when  substituted  into  (F-1)  together  with  T = 2000  K yields 

8 . 

do/o  = 4%,  again  well  below  the  observed  decrease  in  indicated  conductivity. 

The  fact  that  AT  ^ 0.003  T = 6 K justifies  our  assumption  in  the 

o 

previous  paragraph  that  T is  essentially  constant  doim  the  nozzle. 

. ® . . . • • 

Still  another  possible  mechanism  for  lowering  T is  radiation  escape, 

8 

since  the  random  thermal  energy  of  the  gas  must  ultimately  be  the  source  of 
this  energy.  We  may  estimate  the  gas  temperature  drop  which  would  be  caused 
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by  such  a process  by  equating  the  internal  energy  drop  per  second  in  a thin 
slab  of  gas  at  some  location  z to  the  energy  loss  rate  by  radiation  from 
that  slab: 

mde  = R A(z ) dz  , 

Here  m is  the  mass  flow  rate  of  argon,  e is  the  internal  energy  per  unit 

mass,  and  A(z)  is  the  cross-sectional  area  of  the  nozzle  at  axial  position 

z.  Writing  de  = (3/2  h/m  ) dT  ; where  k is  Boltzmann's  constant  and  m 

AS  A 

the  atomic  mass  of  argon,  and  integrating  from  z = 0 to  z = l (the  lengtn 
of  the  nozzle)  yields 


T (0)  - T (L)  = AT 
8 g 8 


R m. 

A 

m ■' 
2 o 


A(z)  dz  . 


By  defining  an  averaged  cross-sectional  area  A such  that 


A L 


fl 


A(z)  dz 


the  previous  result  becomes 

R m A L 

AT  = ^ . (F-11) 

Using  McGregor's  results  as  described  above  to  estimate  R,  and  inserting 
appropriate  estimates  of  A,  L,  m in  (F-11),  we  have  found  that  for  the 
temperature  ranges  of  our  experiments,  AT^  from  equation  (F-11)  is  less 
than  2 K.  The  conductivity  drop  caused  by  this  mechanism  is  therefore  also 
negligible. 
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